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Standard Terminology for Feedback Control 


For a number of years it has been recognized 
that there is a need for standard automatic con- 
trol terminology. In the past, AIEE, ASME, IRE, 
and ASA Committees have been working toward 
this goal in an effort to make technical reading 
and writing easier and to eliminate confusion in 
the meaning of terms which are used to describe 
system performance. 

With the increasing use of feedback control and 
with the demand for higher performance, it is be- 
coming more important that written specifications 
be interpreted properly to ascertain exactly what 
system performance is desired by a customer. 
There have been frequent misunderstandings 
about common terms such as bandwidth, M peak, 
phase margin, damping, rise time, overshoot, set- 
tling time, etc. Actually, a system may be con- 
structed to specifications and yet not meet the 
customer’s performance requirements. Of course, 
it would be desirable to state that the terms re- 
ferred to in a specification are defined in an 
accepted, published standard. 

It is not a function of the PGAC to produce such 
standards for the IRE. This is a responsibility of 
the IRE Committee 26.0 on feedback control and 
its subcommittees 26.1 which considers standard 
terminology for defining a system, and 26.2 which 
considers standards for measuring system char- 
acteristics. Committee 26.2 is relatively new, and 
it has not published any standards for measure- 
ment. However, 26.1 has produced a list of stand- 
ard symbols and terminology which was approved 
by committee 26.0 and eventually by the IRE to 
become the first feedback control standards to 
receive official sanction by a major engineering 
organization. These standards were published in 
the November, 1955 and January, 1956 issues of 
PROCEEDINGS and in the PGAC Transactions of 
February, 1955. They evolved from a concentrated 


study of other existing and proposed standards 
Some were adopted with little change, others were 
modified considerably, and new definitions were 
proposed if necessary. It was the goal of the com- 
mittee to create standards which could be used pre- 
cisely for a large variety of systems, particularly in 
electronic control where the dc loop gain may be 
zero. However, the process of creating standards is 
not easy. It is a difficult, tedious task. Varied opin- 
ions are not easily blended into a satisfactory 
sentence. However, the experience is educational, 
especially when an attempt is made to modify a 
favorite expression of a committee member. 

Although it is difficult for a few committee mem- 
bers to create a definition acceptable to all, it is 
much more difficult to establish a standard which 
will be accepted throughout the control field as the 
definition of a phrase in a textbook or in a specifi- 
cation. For this reason, the committee intends to 
publish future proposed standards on automatic 
control in these TRANSACTIONS solely for the pur- 
pose of provoking comment and constructive criti- 
cism. In this way, it is hoped that the standard 
may be better defined before it is finally approved 
by the IRE feedback control committee and 
adopted as an official IRE standard. The first set 
of proposed standards on terminology is included 
in this issue. Suggestions for rewriting the stand- 
ards, for eliminating them, or for establishing new 
definitions are requested from any interested per- 
son. They should be sent to the editor of these 
TRANSACTIONS where they will be considered for 
publication, if desired, and then forwarded to the 
subcommittee chairman. 

Active assistance in this endeavor is needed, and 
anyone who is interested in establishing new stand- 
ards definitions will be welcome to join subcom- 
mittee 26.1 or 26.2 by contacting the editor of these 
TRANSACTIONS. —The Editor 
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The Issue in Brief 


There are only four papers in this issue of the TRANSACTIONS, 
which is half of the number usually included. It would be interesting 
to know whether readers favor fewer papers in more TRANSACTIONS 
or fewer TRANSACTIONS with more papers. In addition to the papers, 
the annual directory of automatic control engineers is included along 
with Administrative Committee activities and a proposed standard 
terminology for feedback control systems. 

A brief review of each paper follows. 


How the Bandwidth of a Servo Affects Its Saturated Response, George 
A. Biernson—Page 3 


All realizable servomechanisms saturate when large step inputs 
are applied. The response of the servo during saturation differs from 
the predicted linear response, and sometimes this may be detrimental 
to system performance. It is desirable to know when the saturation 
effects may be tolerated. This paper shows not only how the saturated 
response can be predicted from system parameters, but how system 
parameters can be selected to make the saturated response allowable 
for a step input which causes velocity and acceleration saturation. 

By reducing a complex system to a simplified form, the paper 
illustrates that the number of overshoots which result from a velocity 
saturating step response is related to the maximum acceleration, the 
maximum velocity, and the gain-crossover frequency (bandwidth) of 
the primary servo loop. To limit the number of overshoots or to pro- 
duce an optimum response, the bandwidth of the servo can be ad- 
justed. However, if this reduces the bandwidth to an excessively low 
value, it is possible to use nonlinear compensation to allow a stable 
transient for a saturating input and yet maintain a high gain and 
bandwidth for small, nonsaturating inputs. 


Analog Study of Dead-Beat Posicast Control, G. H. Tallman and ~ 


O. J. M. Smith—Page 14 


This paper supplements an earlier paper, “Posicast Control of 
Damped Oscillatory Systems,” by Smith in the September, 1957 issue 
of PRocEEDINGS. Although either paper may be understood without 
reference to the other, the interesting principles of posicast compensa- 
tion can be more easily grasped if both papers are read. The paper in 
this issue describes analog tests on a system where the posicast 
compensation is varied from its nominal theoretical configuration to 
illustrate the effects of imperfect compensation that might be en- 
countered in a practical application. 


Posicast control is also described as a method of splitting the in- 
put excitation into fragments which are delayed in time before being 
applied to the system. When properly used, it eliminates oscillations 
and overshoot in the response of a servomechanism having very 
lightly damped poles. This can be done in a time considerably less 
than one cycle of the uncompensated oscillation. The frequency re- 
sponse of the compensated system can be made flat through and 
beyond the resonant frequency of the system. Compensation is con- 
sidered for simple and complex systems, and for special inputs and 
disturbances anywhere in the forward path of the control system. 


On Closed-Form Expressions for Mean Squares in Discrete-Contin- 
uous Systems, Jack Sklansky—Page 21 


If a system is to be optimized with respect to the mean square of 
some variable, it is usually desirable to obtain a closed-form expres- 
sion for that mean square. This has been an unsolved problem in the 
optimization of discrete-continuous systems where both sampled- 
data and continuous subsystems are used. 

A method of finding the desired closed form expressions is pre- 
sented in this paper. It includes a complete introduction to the sub- 
ject, and the essential steps in the mathematical procedure for analyz- 
ing discrete-continuous systems with random inputs are included. 
The technique is illustrated by deriving closed form expressions for 
the output and ripple of discrete-continuous systems, and for the 
control error of a sampled-data feedback loop. An application to a 
“track-while-scan” system is also included. 


Bibliography of Sampled-Data Control Systems and Z-Transform 
Applications, H. Freeman and O. Lowenschuss—Page 28 


This bibliography represents an attempt to simplify the task of 
searching the literature to obtain a general acquaintance with the 
field of sampled-data control and to obtain information on certain 
specific aspects of it. Due to the close association of the Z-transform 
and sampled-data systems, a number of papers dealing exclusively 
with the Z transform are included in the bibliography. 

The bibliography includes only that material which in the au- 
thors’ opinion represents either a significant contribution to the field 
or has tutorial value. The bibliography is arranged alphabetically 
according to the name of the author, or the first-named author in 
cases of coauthorship. A subject index with cross references to the 
author list is also provided. 
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How the Bandwidth of a Servo Affects 


Its Saturated Response’ 
GEORGE A. BIERNSON+ 


Summary—The transient response of a servomechanism to a 
very large step may be quite oscillatory because of saturation of the 
power member in acceleration. By making certain idealized yet 
reasonable assumptions, this paper shows that the number of over- 
shoots of the response to a step large enough to drive the system 
to velocity saturation is related quite simply to the maximum ac- 
celeration, maximum velocity, and the gain crossover frequency 
(bandwidth) of the primary servoloop. To limit the number of over- 
shoots in certain applications, the maximum velocity and/or band- 
width of the servo may have to be excessively lowered. However, it 
is possible to remedy this difficulty with nonlinear compensation, 
which varies the gain crossover frequency inversely with the output 
velocity, so as to allow a stable transient for a large step input yet 
maintain high gain and a wide bandwidth for small inputs. 


INTRODUCTION 


N important nonlinear problem encountered in 
many feedback control systems is a poorly 
damped oscillating transient following a very 

large step input, which is due to saturation of the second 
derivative, 7.e., saturation of the acceleration of the 
controlled variable. This saturated step response occurs 
in many servomechanisms when the system is syn- 
chronizing on a particular reference input signal. 

Although the saturated step response may be non- 
linear, the important linear parameters of the system 
are nevertheless of fundamental importance in deter- 
mining its form. The purpose of this paper is to relate 
in a simple fashion the form of the saturated step re- 
sponse to the basic linear parameters of the system and 
to the saturation levels of the controlled variable. 

The most important linear parameter of the system 
is the gain crossover frequency, defined as the frequency 
at which the magnitude of the open loop transfer func- 
tion, or loop gain, is unity. For the linear step response, 
the time to rise to within 37 per cent of the final value 
is roughly equal to the reciprocal of the gain crossover 
_ frequency. Since the loop gain of most feedback control 

loops approaches an asymptote inversely proportional 
to frequency in the region near gain crossover, it often is 
desirable to consider, instead of the gain crossover fre- 
quency, the more convenient parameter, the asymptote 
crossover frequency. That is, the frequency at which 
that asymptote is unity. Generally, the asymptote 
crossover frequency is closely equal to the actual gain 
crossover frequency, and consequently the asymptote 
crossover frequency is employed here, but is referred to 
for convenience as the gain crossover frequency and 
designated as w,. 

To make the concepts presented here more concrete, 


* Manuscript received by the PGAC, May 5, 1957. 
+ Sylvania Electric Products, Inc., Waltham, Mass. 


the analysis is developed in terms of a positional servo- 
mechanism, but nevertheless can apply to any feedback 
system. The systems studied have saturation in both 
velocity and acceleration, which in a general sense repre- 
sents saturation in the first and second derivatives of 
the controlled variable. 

The impetus for this investigation was a linear anal- 
ysis of the synchronizing transient of servomechanisms 
made by Travers.! He showed that in order for the 
drive members of a servo to keep from saturating in ac- 
celeration at the end of a large synchronizing transient, 
the maximum value of the system impulse response 
should be less than the ratio Am/Vm, where Am is the 
maximum motor acceleration and V,, the maximum 
velocity. Since for well-designed systems the maximum 
value of the impulse response is roughly equal to the 
gain crossover frequency? w-, this gives roughly the fol- 
lowing requirement for linear deceleration of the syn- 
chronizing transient. 

Gea an . Gi 

Vian 

Of course, it generally is desirable that some accelera- 
tion saturation occur at the end of the synchronizing 
transient in order to use the drive motor effectively. 
However, Travers reasoned intuitively that if the gain 
crossover frequency were increased more than a factor 
of ten greater than the ratio Am/Vm, a very poor syn- 
chronizing transient would probably result. The anal- 
ysis of this paper shows that Travers’ reasoning was 
essentially correct. 

The material of this paper was presented originally 
in a report by the author.’ 


ANALYSIS OF A RATE FEEDBACK SYSTEM 
Description of System 


The first system to be considered is the rate feedback 
system shown in Fig. 1. It is a good approximation to a 
large class of feedback control systems. It could repre- 
sent a positional servomechanism with tachometer 
feedback. The variables and parameters of the system 
are listed. 


1P, Travers, “Motor Saturation in Servomechanisms,” Mass. 
Inst. Tech., Cambridge, Mass., Dept. of Elec. Eng., Master’s Thesis; 
September, 1948. 

2 G. A. Biernson, “Estimating transient responses from open-loop 
frequency response,” AIEE Trans., Applications and Industry, pp. 
388-403; January, 1950. 

3—, “Relating the Saturated Step Response of a Servo to Its 
Bandwidth,” M.I.T., Cambridge, Mass., Servomechanisms Lab., 
Rep. 7138-TM-2; June 3, 1955, 


Fig. 1—Rate feedback system. 


X,=reference input (position). 

X,.=controlled variable (position). 

X,=error in position. 

V.=velocity of controlled variable. 
A,=acceleration of controlled variable. 
@,=gain crossover frequency of position loop. 
Wer = gain crossover frequency of rate loop. 


The absolute values of speed and acceleration are 
designated by the symbols V and A: 


V=|Vel, (2) 

A=|A,|. (3) 

The maximum speed and acceleration capabilities of the 
controlled variable are designated as Vm and Am: 
Vn = maximum (V) = maximum | V,|, 
Am = maximum (A) = maximum | A. | 4 

The significance of the frequencies w, and we, is il- 

lustrated by the plots in Fig. 2. The open loop transfer 


function of the rate loop (with the position loop opened) 
is shown by Fig. 1 to be 


Go) = Grol 5: (4) 
The closed loop response of the rate loop is 
Gor) Wer 
OK : (5) 
A Gor) S =P Wer 


The open loop transfer function of the position loop 
(with the rate loop included) is 

We Gor) We 
lata © 
s LI1+ Gey s[(s/@er) + 1] 
Frequency response plots of | Gon| and | Gos are shown 


in Fig. 2. These are plots of the loop gains of the posi- 
tion and rate loops. 


Gp) a 


Ideal Saturated Step Response 


The ideal saturated step response can be considered 
one in which the controlled variable is saturated at all 
times either in acceleration or velocity and which has 
no overshoot. Such a transient yields the fastest pos- 
sible synchronizing response for the given saturation 
limits of the controlled member. The gain crossover fre- 
quencies w, and w,, required to give this ideal transient 
shall be determined now. 

Fig. 3 shows both the position X, (or error X,) and 
velocity V. for the ideal nonlinear response of the sys- 
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Fig. 2—Loop gains of position and rate loops. 


REFERENCE 
INPUT X, 


POSITION X.—7 > 


VELOCITY V, 


TIME 


Fig. 3—Ideal saturated step response. 


tem of Fig. 1 toa large step input. At time 4; the step is 
applied, and between #4 and ¢ the controlled member is 
under maximum acceleration A,. At ft. the controlled 
variable reaches maximum velocity Vi», and between fp 
and ¢; it moves with constant velocity. Now in order for 
the transient to end in the shortest possible time, it is 
necessary for the system to start decelerating at a point 
3, such that the velocity is reduced to zero at the time 
the error reaches zero, at point 4. Thus, between points 
3 and 4 the controlled member should be under a satu- 
rated deceleration of Am. 

The time for deceleration (t,;—#;) obviously is equal 
to the time for acceleration (t:—4), which is equal to 


(i=) = CG eis i | 
rot Foal eta _ (7) | 
The motion of the controlled variable during accelera- 
tion and during deceleration is equal to (1/2)An?, ; 
which is equal [by (7) | to : 
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Vin ; 
: 


aD (8) 


1 1 aaNe 1 
— Ant” An( ) = 
2 2 Am 2 
This distance is shown in Fig. 3. Eq. (8) thus shows that, 
in order to achieve the ideal saturated step response for 
an input large enough to produce velocity saturation, 
the system must “switch” discontinuously to maximum 
deceleration at an error X, equal to Vn?/2Am. 

The manner in which the system switches to maxi- 
mum deceleration can be seen by examining Fig. 1, 
which shows the acceleration to be equal to 


A = wer(weXe — Ve). (9) 
Eq. (9) shows that the acceleration is positive when 
aekee> Ve (10) 
and it is negative when 
He eR (11) 


Thus, the system switches between positive and nega- 
tive acceleration at an error of 


X. = Ve/ee. (12) 


If w., is made much greater than (Am/Vm), (9) shows 
that a very small change of X, from the value of (12) 
produces either maximum acceleration or deceleration 
(i.e... A=+A, or A=—A,). Thus, for large values of 
the rate loop crossover frequency a ,, the system 
switches almost instantaneously between maximum ac- 
celeration and maximum deceleration at the value of 
error given by (12). 

Prior to the point of switching to deceleration, the 
controlled variable speed is equal to Vn, so that the 
switch to maximum deceleration occurs at the error 
Vm/wc. Therefore, in order that the saturated step re- 
sponse not overshoot, this switching must take place at 
an error larger than the distance required for the con- 
trolled member to decelerate to zero speed, which is 
Vn? /2Am. Hence for no overshoot 


Ve Wit 
(13) 
We DAs 
which gives 
245 
we <<—— - (14) 


For the ideal saturated step response, w, should be equal 


2A el Vines 


Saturated Step Response with Overshoot 


If w, is increased above the value 2(Am/Vm), over- 
shoot occurs in the saturated step response. An analysis 


-presented in Appendix I shows that the number of over- 
shoots x of the saturated step response is equal to. 


e's, 


Vite (15) 
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When 1 is not an integer, the number of overshoots is 
the next highest integer. Solving (15) for w, gives the 
convenient expression 


7 A m 


Op ar) (16) 


When 1 is set equal to zero, (16) reduces to the optimum 
value for w,., which is 24 m/ Vm. 
The settling time of the saturated step response is 


LT, = Vm/Am)(1 + 1). (17) 


This does not include the time during which the con- 
trolled member is moving at saturated velocity (the 
time between points 2 and 3 in Fig. 3). The minimum 
size of the step input required to produce a full saturated 
step response (7.e., to produce saturated velocity) is 


Ving 
2Am 


+2 
Xoo (eet oy ¢ ) (18) 


n= 1 
The peak overshoot is related to this minimum input 
step by 


Pow = [n/(n + 2)|Xereat- 


The above equations apply strictly only for integer 
values of x. A noninteger value for 2 indicates that in 
the last segment of the transient the system operates 
linearly with respect to the position loop, 7.e., with a 
velocity proportional to error rather than with a con- 
stant acceleration. 


(19) 


Effect of Rate Loop Crossover Frequency 


The above analysis has assumed that the gain cross- 
over frequency of the rate loop w.; is very much greater 
than the ratio (Am/Vm): 


wy > (Ag/ Va): 


Substituting (16) into (20) gives an alternate expression 
for this requirement: 


@epi@y/ 21-1): (21) 


An approximate analysis of the effect of a low value of 
Wer is made in Appendix II. The analysis shows that the 
value of w, for no overshoot is approximately 


pies) 


If w., is infinite, w, is equal to the value 24 m/Vm calcu- 
lated previously in (16) when z is set equal to zero. As 
Wer is decreased from infinity, w. must be decreased if 
there is to be no overshoot. 

Usually the gain crossover frequency of the rate loop 
is at least twice that of the position loop, 


(20) 


We 


(22) 


du cr 


Wer = 2a. (23) 


For this condition, the value for w, in (22) is no less than 
7/8 of the value it has for infinite w.,. Thus, as long as 
the rate loop has a gain crossover frequency at least 
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twice that of the position loop, one can assume that 
Wr is infinite without making much error in calculating 
the characteristics of the saturated step response. 


Effect of Other Dynamic Elements in the System 


Additional high-frequency lags in the systems, which 
always occur in practice, would not alter the saturated 
step response appreciably, as long as the systems are 
adequately stable under linear conditions. The primary 
effect of such lags on the saturated step response is to de- 
lay the switching and hence somewhat prolong the 
transient. 

On the other hand, low-frequency integral networks 
in the system have a very pronounced effect upon-the 
saturated step response. This effect in some cases is so 
bad that it is necessary to short the integral networks 
while synchronizing. Such shorting, of course, does not 
limit the accuracy of the system, because integral net- 
works are not designed to operate under such conditions. 
Leaving the integral networks in the system during this 
transient not only tends to produce instability, but also 
adds a tail to the transient. 


EXTENSION TO OTHER SYSTEMS 


Besides the rate feedback system described in the 
preceding section. there are two other common types of 
systems to which the analysis can apply if suitably 
modified. These are the velocity source system illus- 
trated in Fig. 4 and the lead network system in Fig. 
5(a). The lead network system uses a lead network for 
compensation rather than rate feedback. The velocity 
source system needs no compensation, because the 
power member provides a velocity proportional to 
applied signal. An example of a typical velocity source 
system is a hydraulic servo employing a valve which 
gives an output flow proportional to valve displacement. 


Velocity Source System 


For the system of Fig. 4, the nondimensionalized 
variable E,, represents the signal applied to the motor. 
The torque-speed curve of the motor shows that for a 
given value of £,, the motor supplies whatever accelera- 
tion (z.e., torque) is required, up to its maximum accel- 
eration capability, to drive the controlled variable at a 
velocity proportional to E,,. Therefore, the system acts 
just like a rate feedback system (Fig. 1) with an infinite 
gain crossover frequency in the rate loop. The velocity 
source system is under maximum positive acceleration 
whenever w,X. is greater than V,, and is under maxi- 
mum negative acceleration when w,X, is less than Vp. 
Hence, the equations that have been derived for the 
rate feedback system of Fig. 1 apply directly to the 
velocity source system of Fig. 4. 


Lead Network System 


In Fig. 5(a), the lead network transfer function is 
shown as (s+w ,), which approximates the actual trans- 
fer function of a lead network: 
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Fig. 4—Velocity source system. 
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Fig. 5—Lead network compensated systems. (a) Lead network 
system. (b) Equivalent rate feedback system. 


stor (/a)(s+ a) 


(s/awy =in 1) (24) 


S+ aw) 


in cascade with a gain of a. The block diagram neglects 
the factor [(s/aw:)+1], but this factor has little effect 
upon the saturated step response if a is reasonably 
large. 

An analysis of the lead network system is presented in 
Appendix III. The analysis shows that during the 
saturated step response the lead network system of 
Fig. 5(a) behaves essentially the same as the equivalent 
rate feedback system in Fig. 5(b). The equivalent rate 
feedback system has a position loop crossover fre- 
quency equal to w, and a rate loop crossover frequency 
equal to w;. Applying (16) and (23) derived for the rate 
feedback system gives 


be > Dan. (26) 


; 


q 
4 
' 


/ 


L958. 


The condition of (26) is almost always satisfied in prac- 
tice. 

Thus, the equations derived for the rate feedback 
system apply to the lead network system when w, for 
the rate feedback system is replaced by the lead net- 
work break frequency a7. 


NONLINEAR TECHNIQUES FOR 
IMPROVING PERFORMANCE 


It has been shown that limitations must be placed 
upon the position loop gain in order to achieve adequate 
synchronizing performance. In many cases, however, 
these restrictions result in a great decrease in the linear 
performance of the system, and hence it is very desir- 
able to employ nonlinear techniques to allow both a fast 
linear response and an adequately damped saturated 
step response. 

The principle for achieving this is based upon the 
equation for the ideal saturated step response obtained 
by setting 2 equal to zero in (16). 


2Am 
aN 


2, = (27) 
The equation shows that, for an ideal step response to a 
saturating step input, the gain crossover frequency 
should be inversely proportional to the maximum speed 
of the controlled variable. Hence, ideal performance 
can be obtained for all operating speeds if w, is varied 
inversely with the speed. Of course, this variation of w, 
could not be achieved for very small speeds because 
it would require an excessively large w., which would 
lead to instability because of linear effects. Conse- 
quently, when the system is operating at low speed, w, 
should be set at the maximum value commensurate 
with adequate linear stability, and when the speed in- 
creases above the value 2A »/w, the frequency w, should 
be decreased inversely proportional to the speed. 

A simple way to instrument this change of w, in a 
rate feedback system is to place a quadratic gain ele- 
ment in the rate feedback path,‘ so that the gain in this 
path is proportional to the magnitude of the velocity 
signal. This has the effect of changing the gain in the 

‘loop (and hence the gain crossover frequency w,) in- 
versely with the controlled variable speed. On the other 
hand, it has the disadvantage of varying the gain cross 
over frequency in the tachometer loop, so that, if the 
tachometer loop is to be kept stable at high speed, it 
must operate at excessively low bandwidth at low speed. 

Another way of improving the stability of the satu- 
rated step response without excessively lowering the 
gain crossover frequency is to limit the maximum motor 
speed V,, by placing a governor on the motor.® This, of 
course, has the disadvantages of greatly slowing down 


_ 4W.H. Surber, Jr., “A nonlinear compensating configuration for 
saturating servomechanisms (abstract),” 1955 IRE CONVENTION 
-Recorp, Part 4, p. 23. at ; 
8 W.C. Robinette, “A packaged servomechanism,” Electronics, 


" 1948. 
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the synchronizing transient, but it affords a simple 
means of stabilizing a low-performance feedback con- 
trol system. 


EFFECT OF MOTOR BREAK FREQUENCY ON 
SATURATED STEP RESPONSE 


The equation relating the gain crossover frequency to 
the number of overshoots on the saturated step re- 
sponse can be expressed in terms of the motor break 
frequency (or its reciprocal, the motor time constant), 
provided the motor break frequency (or time constant) 
is defined on a nonlinear basis. To show this, it is nec- 
essary to review the commonly used linear definition 
for the motor break frequency. 

If the torque-speed curve of the motor is linear as 
shown in Fig. 6(a), (28) holds for a constant value of 
motor voltage E,. 


jpn PEE EU 
=e Ga iE EAS 


where 7, is the stalled torque (torque at zero speed), 

Jm is the motor inertia (including load), and f, is the 

negative of the slope of the torque-speed curve, that is, 
Ga 


fm = — — for constant E,. 
dV 


(28) 
(29) 


(30) 


The stalled torque generally is nearly proportional to 
the applied motor voltage E,,, so that one can assume 
the relation 


le — TEs (3 1) 


where the proportionality constant K,, may be called 
the motor gain. Combining (29) and (31) gives for the 
transfer function between applied voltage and motor 
speed 


V En 
Aish eS a Cee (32) 
En 2 Sm Aes 
By defining the motor time constant T,, as 
Tm = wikia (33) 


this transfer function may be expressed in the non- 
dimensional form 
Vo Kalla | ay 
VR get et | 
In many cases, however, the torque-speed curves are 
not linear and hence the above analysis does not strictly 
apply. For example, Fig. 6(b) shows a typical torque- 
speed curve of a two-phase servomotor driven from a 
high-impedance source, such as a pair of 6L6 vacuum 
tubes. In considering the transfer function of the motor 
one must generally include the power amplifier and 
motor together as a single unit, because of interaction 
between the two, and hence the voltage E,, should be 


taken as the input signal to the power amplifier rather 
than the voltage applied to the motor. The torque-speed 
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Fig. 6—Torque-speed curves. (a) Linear torque-speed curve. (b) Nonlinear torque-speed curve. (c) Definition of saturated damping factor. 
d) Comparison of decelerating torque with stalled torque. 


curves of a two-phase servomotor driven from a low- 
impedance source would be much more linear than 
those of Fig. 6(b). 

For nonlinear torque-speed curves the same transfer 
function derived on a linear basis is still used, but the 
damping ratio fm, defined as the slope of the torque- 
speed curve, is considered as a variable depending upon 
where the motor is operating. Since operation near zero 
speed is often considered to be the most important, the 
damping factor generally is taken as the negative of the 
slope of the torque-speed curve at zero speed. 

On the other hand, this linearizing approach leads 
to the following paradox. The gain factor K,,/fm in (34) 
is inversely proportional to the motor damping and 
hence is inversely proportional to the negative of the 
slope of the torque-speed curve. For the torque-speed 
curve in Fig. 6(b), this slope can vary from zero at low 
speed to a fairly large negative value at high speed, and 
hence the gain factor K»,/fm of the motor can vary from 
a rather low value at high speed to infinity at low speed. 
Since this gain factor is a cascade element in the loop, 
this great variation in motor damping must produce a 
great variation in the loop gain, which may lead one to 
wonder how the system can have reasonable dynamic 
performance for all values of motor speed. 

The answer to this paradox can be seen quite simply 
by rewriting the motor transfer function of (34) in the 
form 


V Ki Sin 
S + wm 


(35) 


where w,, is the motor break frequency, which is the 
reciprocal of the motor time constant Tm: 


1 Tin 
Cy are 5 Seek 
Tm Te 


(36) 


The most important transient test of the dynamic per- 
formance of a feedback control loop is its step response, 
and the rise time and peak overshoot (or stability) of 
the step response are determined primarily by the char- 
acteristics of the loop transfer function in the region 
near gain crossover. Hence variations in the motor 
transfer function V/E,, are of basic importance only if 
they occur in the frequency region near gain crossover. 
Systems with torque-speed curves like Fig. 6(b), gen- 
erally have compensation to provide system damping, 
so that the resultant gain crossover frequency w, is much 
larger than the motor break frequency w,. For this con- 
dition the motor transfer function of (35) can be closely 

approximated as 
V Kn/Im — Km/Im — Km/Im 


Em jo + om 
in the frequency region near gain crossover. Thus, varia- 
tions in the motor damping are of little importance, be- 
cause the motor break frequency has negligible effect 
on the loop transfer function near gain crossover. 

Of course, if the system does not have additional com- 
pensation to give damping, the motor damping is very 
important. For such cases the approximation of (37) 
cannot apply at gain crossover, because instability re- 
sults if the gain crossover frequency is made much 
larger than the motor break frequency wm, and there- 
fore it often is more convenient to express the motor 
transfer function in terms of the motor time constant 
Tm as in (34). On the other hand, when there is no damp- 
ing supplied by compensation, the motor damping must 


(37) 
Jo 
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be quite large for all values of speed and rather con- 
stant in order to achieve reasonable dynamic perform- 
ance. For this condition, large nonlinearity of the motor 
damping generally is not encountered in practice. 

For wide bandwidth systems in which the damping is 
provided by compensation, a variation in motor damp- 
ing merely changes the velocity constant, that is, it 
changes the steady-state errors due to constant velocity 
inputs. For the purpose of minimizing these steady- 
state errors it is desirable that the motor damping be as 
low as possible. If the motor damping is zero at low 
speeds, the velocity constant is infinite, and the steady- 
state errors for low velocity inputs are zero. 

Thus, it has been shown that the term “motor-break- 
frequency Wm,” or its reciprocal, the “motor time-con- 
stant Tm,” defined on a linear basis, has little effect upon 
system performance unless the system is operated at a 
rather low bandwidth, with the motor damping furnish- 
ing the main damping for the system. On the other 
hand, there is a similar parameter defined on a nonlinear 
basis which is quite important, and this parameter shall 
be called the saturated motor break-frequency Wms. 

The definition for the saturated motor break-fre- 
quency wn, is illustrated in Fig. 6(c); it shows the torque- 
speed curve for the maximum voltage E,,. The maximum 
motor speed is shown as V,, and the maximum motor 
torque as 7;,,. Ifa straight line is drawn from the points 
of maximum torque and maximum speed, as shown, the 
negative of the slope of this line is defined as the satu- 
rated motor damping factor fms, so that 


Tims re AE 


The saturated motor break-frequency Wms then may be 
defined as 


(38) 


Toma sie Be 
eee 2 


Wms (39) 

One can also define a saturated motor time-constant 
Tms as being equal to the reciprocal of wns. Newton® has 
defined the “motor time-constant” in this manner, and 
has shown its importance in describing the adequacy of 
the power capabilities of a motor in various applications. 
He points out that the saturated motor time-constant 
Tms Tepresents the time for the motor to accelerate from 
zero to maximum speed. 

By means of the basic relation of (16) between w, and 
the number of overshoots of the saturated step response, 
the gain crossover frequency w, can be related to the 
saturated motor break frequency. For positive speeds, 
the maximum motor acceleration occurs at zero speed 
and is equal to 


_ However, in terms of the saturated step response it is the 
deceleration capabilities of the motor that are of prime 
importance, that is, the acceleration at negative speeds. 


6 G. C. Newton, Jr., “What size motor?” Machine Design, No- 
at Behe, 1950. 
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If the torque-speed curve is as shown by curve 1 of 
Fig. 6(d), the maximum deceleration is the same as the 
value given in (40), but if it is as shown by curve 2 the 
maximum deceleration is somewhat larger. However, 
for simplicity, one can take (40) as representing a rea- 
sonable approximation to 4,,. Combining (39) and (40) 
relates the maximum acceleration to the saturated 


motor break-frequency by 
Onin Vex (41) 


Now, the basic equations for the saturated step response 
have been given in (16) and (25) as 


2Am rate feetback system 
an == (1+ n) : (42) 
an velocity source system, 
DAT 
a, = 7 (1 + n) lead network system. (43) 
Substituting (41) into the above gives 
rate feedback system 
(Oya 2(1 SF 1) Wms (44) 
velocity source system, 
w, = 2(1 + 2)wme— lead network system. (45) 


These equations thus show that the saturated motor 
break frequency is an important limitation on system 
bandwidth. 

It is interesting to point out that if the torque-speed 
curves are linear as shown in Fig. 6(a), the saturated 
motor break-frequency @ms is equal to the linear motor 
break-frequency w,. This illustrates why the motor 
break-frequency w, has been found to be such a limita- 
tion on system bandwidth when the torque-speed curves 
are linear, even when additional compensation is added 
to provide damping. 


CONCLUSION 


Because of saturation of the second derivative (ac- 
celeration) of the controlled variable of a system, the 
response to a very large step input can be very poorly 
damped, and may even be unstable. This nonlinear in- 
stability is due to the inability of the motor to deceler- 
ate the controlled variable adequately as the error 
reaches zero, so that the controlled variable is driven 
well into saturation in the reverse direction before 
coming to a stop. To remedy this problem, either the 
maximum speed of the controlled variable must be 
limited or the loop gain must be decreased, in order that 
the controlled variable start decelerating at a larger 
value of error. 

In systems where instability of the saturated step re- 
sponse is a problem, nonlinear compensation can be em- 
ployed to allow both high gain during linear operation 
and a fast well-damped response to large saturating step 
inputs. This is done by varying the loop gain inversely 
with the speed of the controlled variable. 

For a system employing a velocity-source power 
member (such as is used in many hydraulic servomech- 
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anisms) the gain crossover frequency w, is related as 
follows to the number of overshoots of the step re- 
sponse for an input large enough to produce maximum 
saturation: 


Dales 


46 
_ ), (46) 


W. = 


where A,, is the maximum motor acceleration and V,, 
the maximum velocity. This same relation also applies 
to rate feedback compensated systems, employing an 
acceleration (torque) source power member, if the rate 
loop gain crossover frequency wr is much larger than 
the position loop gain crossover frequency w,. Excluding 
the time the controlled variable is slewing at saturated 
velocity, the setting time for the transient is 


2 Vin 
Rs Day BACs ast i era (47) 


We m 


For lead network compensated systems, in which the 
gain crossover frequency is significantly greater than 
the lead network break-frequency w;, the above rela- 
tions apply if w, is replaced by a. 


APPENDIX I 
ANALYSIS OF OSCILLATORY STEP RESPONSE 


When the gain crossover frequency w, of the rate 

feedback system is made much greater than 24m/Vm, 
the saturated step response becomes quite oscillatory. 
This Appendix presents an analysis of this oscillatory 
step response. 
» Fig. 7 shows an oscillatory saturated step response of 
the rate feedback system of Fig. 1. Points 2, 4, 6, and 8 
represent the points where switching of acceleration 
occurs, and hence represent points of maximum speed; 
while points 1, 3, 5, 7, and 10 represent points of zero 
speed. For this transient the value of w, has been set so 
that following the last point of acceleration (point 2) the 
system decelerates with maximum deceleration, reach- 
ing zero velocity at zero error. 

Jn the analysis to be presented, the following symbols 
are employed: 


X,=value of X, at point x. 
Xen) =Value of X, at point x. 
V,=speed at point n. 


At points 2, 4, 6, and 8, where the acceleration switch- 
ing occurs, (12) shows that the magnitude of error is re- 
lated to the speed by 


| Xeim| = Vn/ eve. (48) 
This gives 

| Xe) | = V2/We (49a) 

|X oc) | = Vales (49b) 

| Xe) | = V6/w (49c) 

| Xe) | = Ve/we. (494) 


March 
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Fig. 7—Saturated step response with overshoot: 


Fig. 7 shows that the magnitudes of error at successive 
switching points are related as follows. 


| Xea| = | X2— Xi (50a) 
| Xea| = |Xe—Xs| — | Xs— X2| — | Xew| (SOb) 
| Xe)| = |Xs— Xs| — | X¥s— X4| — | Xe] (50) 
| Xe) | = |Xs— X7| — | X7 — Xe| — | Xece)|. (50d) 


During the entire transient between point 8 and 
point 1, full acceleration of magnitude Am is applied, 
alternately positively and negatively. Hence the dis- 
tances AX between successive points of zero and maxi- 
mum speed are equal to 
1 ave 
— A,,(At)? = — —, 
u ra? 


™m 


AX = (51) 
where V is the maximum speed for that section of the 
transient. Applying (51) to the specific regions of the 
transient gives 


Li Vat 

Pe ee hone eae (52a) 
2A 
gor 
2 Aas 
2 ¥3° 

| ety | corel Mp eG lige nea he 
2) Age 
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|Xs— X7| = |Xw— Xs 


=——.- (52d 
To (52d) 


Substituting (50a)—(50d) and (52a)—(52d) into (49a)— 
(49d) gives 


Ve Ue qice 
= (53a) 
We Li. Ae 
Va te rg PP ie V2 
= — = (53b) 
We 2 ae 2 An We 
Ve i Ve i V 22 V4 
=— — — — — — (53c) 
We 2 A 2 Aw We 
Vs 1 V3" i! Ve" Ve ee 
= — _ _ . (53d) 
We Daa AS, Die As We 
Rearranging the individual equations gives 
2Am 
Vo. = (54a) 
We 
ZA se eA. 
VA2 ae ie | Vo — ve | = 0 (54b) 
We We 
DA en 2Ax 
Ve, 7a 65 | Vs ae ve = 0 (54c) 
De We 
DA ws 2Am 
Ve? — Vs — | Ve + Vo? | = 0. _(S4d) 
We We 


Solve (54b) for V4, 


rote (Ere) 


The polynomial under the square-root sign is a perfect 
_ square. Choose the plus sign, because the speed V4 must 
be positive. This gives 


Aen Aim 2Am 
Va= + —+ V, = —+ D2. (56) 
We We We 
Thus, the solution of (54a)—(54d) is 
Am 
Vo = 2 — (57a) 
We 
Vi = +V2= 2( ) (57b) 
We We 
DAG QA 
eV ge 3( ) (57c) 
We We 
DAy DAs 
We We 


Now the speed Vs is the maximum speed of the sys- 
tem V,,. Hence (57d) gives for this transient, which has 
_ three overshoots, 

/ 2Am 
Vin = 


& We 


(4). (58) 


Hence, for the case of 2 complete overshoots, 


Zain } 
Van = (wm + 1). (59) 
We 
Solving for w, gives 
Get 60) 
On n ; 
Vn \ 


Eq. (60) is thus a generalization of (14), which applied 
to the ideal saturated step response. Eq. (60) shows that 
the greater w, is above 2A m/Vm, the greater is the num- 
ber of overshoots in the saturated step response. Solv- 
ing (60) for 2 gives 


VinWe 
a ip 


n Ls (61) 
This equation, of course, applies strictly only for integer 
values of x. A noninteger value for indicates that in 
the last segment of the transient the system operates 
linearly with respect to the position loop, that is, with a 
velocity proportional to error rather than with a con- 
stant acceleration. 

The settling time of the transient can be computed as 
follows. Designate the time interval between two suc- 
cessive points a and 0 as fq. The time intervals of the 
various segments of the transient are equal to 


tio = tos = Vo/Am (62a) 
tsa = tas = Va/Am (62b) 
tss = ter = Vo/Am (62c) 
trg = to-10 = Vs/Am. (62d) 


Neglecting the time the controlled variable is slewing 
at maximum velocity, between points 8 and 9, the total 
settling time is 


es = ie + oe ye (63) 


2 

Sa ay hl ont) (64) 
4 

x (ote 32" a)! (65) 
We 


For a transient with nm overshoots, the settling time is 
thus equal to 
4 


T,.=—[1+2+-+---+n4+ (n+). (66) 


We 


Now the sum of integers up to q is equal to 


= 1) 
142434 ---ga Me’, (67) 
Hence (66) is equal to 
y} 
T, = — (n+ 1)(n + 2). (68) 


We 
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In terms of the saturation parameters, the settling time 
is equal to 


2) Van 
Toe ae (69) 


It is also important to know the size of the overshoot 
in the transient. The magnitudes of error at the points 
of zero velocity are equal to 


eee 70 
| Xeca | ny Gaae (70a) 
1a Vacamen yc 
35 70b 
SG Serene . (70b) 
alee ate 
> ete 70 
| Xe | ae a . (70c) 
tae 
og |, X,|. (70d) 
A We 


Substituting (57a)—(57d) for the speeds in the above 
gives 


DvAn 

PX eae ms — (1b) (1) (71a) 
(Shoo 
DAT, 

| Xe = —— (1 + 2)(2) (71b) 
We 
2, 

| Xem| =— A + 3)@) (71c) 
We 
2Am 

| Xecao | (1+ 4)(4)+ AX, (71d) 


where AX,,4 is the motion of the controlled variable 


while under saturated velocity, given by 
XM re, — Xa | (72) 


Eq. (7ic) represents the peak overshoot Py. For the 
general transient of 2 overshoots, the peak overshoot is 
thus 


2A 


Poy = —— n(1 +n). (73) 

We" 

More convenient expressions are 

NV m, 
Pee I (74) 

We 
( s ) 75 
ao Gott) | 75) 


Eq. 74 shows that the peak overshoot is 7 times the 
maximum linear range of the position loop, which is 
Vin/we. Eq. (75) shows that the peak overshoot varies 
with 2 by the factor n/(w+1). Thus, as long as the 
transient has at least a couple of overshoots, the peak 
overshoot does not increase much with increasing n. 


March 


The maximum error occurs at the start of the tran- 
sient [representing X-ao) for the example | and is equal to 


DeAan 
(n + 1)(m + 2) + AXzat. (76) 


2 
We 


X eGane) = 


If AXsat, the displacement during maximum velocity, is 
zero, (76) gives the error for the smallest step which will 
bring the controlled variable up to saturated velocity, 
which is 


LAs 


2 


(1 + 2) Vm? 


(pi) 2A (77) 


(1+ 1)(w+ 2) = 
WO 

This varies only by a factor of two as the number of 

overshoots 7 is increased from zero to infinity. 


APPENDIX II 


EFFECT OF GAIN CROSSOVER 
FREQUENCY OF RATE LOOP 


The analysis of the rate feedback system presented in 
Appendix I assumed that the gain crossover frequency 
of the rate loop was infinite. This Appendix develops a 
simple approximation of the effect of lowering the rate 
loop gain crossover frequency. 


Fig. 1 showed that the-acceleration A, is equal to 
Ap. = WelerXo— Ve). (78) 


Hence, the controlled member is under maximum posi- 


‘tive acceleration when 


Am 


Ve 
YES Brae 


We Weer 


(79) 


and it is under maximum negative acceleration when 


Ve a 


We WeWer 


Xe (80) 
Fig. 8 shows approximately the saturated step re- 
sponse which occurs with a lower w-,, if w, is adjusted for 
the maximum value with no overshoot. At the error 
Vin/@e (point 3) the acceleration begins to become 
negative, but it is not until point 5 that the full negative 
acceleration (—A,) is applied. Hence, at point 5 the 
error is 
Vs Am 
AY elie ae : 


We WeWer 


(81) 


If it is assumed that there is negligible change in velocity 
between points 3 and 5, the error at point 5 is approxi- 
mately equal to 


(82) 


Thus, lowering the rate loop gain crossover frequency 
decreases the error at which maximum deceleration is 
applied by the amount Am/wWer. 

If the velocity is assumed nearly constant between 
points 3 and 5, the acceleration increases essentially 


: 
; 
| 


J 


1958 Biernson: How the Bandwidth of a 


INPUT X, 


POSITION X, 


ACCELERATION A, 


Fig. 8—Effect of rate loop gain. 


linearly with the displacement in this region, as well as 
with time. The approximate acceleration plot is shown 
in Fig. 8(b). The change in velocity between points 3 
and 5 is equal to the area under the acceleration curve 
in this region, which is equal to 
1 
AV = ip A dt = — Am(ts — 13). (83) 
ts y, 
This same change in velocity would occur in this region 
if the acceleration were zero up to point 4, and were 
equal to A,, thereafter, where point 4 is midway between 
points 3 and 5. This acceleration curve is shown by the 
dashed segment in Fig. 8(b). 

Thus, a reasonable approximation of the transient can 
be obtained by assuming that the acceleration switches 
discontinuously from zero to maximum deceleration at 
point 4. Since point 4 is midway between points 3 and 5, 
this effective switch of acceleration occurs at the error 


(84) 


We 


If the controlled variable is to decelerate from point 4 to 
zero velocity at zero error, the error at point 4 must be 
equal to 


ee (85) 
Xow = — 8 
Nee A. 
‘| Equating the expressions of (84) and (85) gives 
2Amn 2Am/ Vin 
w= {1 - “|, (86) 


Aer 
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Since (24n/Vm) is roughly equal to w,, (86) is roughly 


equal to 
(et 


Eq. (87) shows that, as a first approximation, the 
effect of decreasing w,, is to require that w, be decreased 
by the factor [1—(w,/4w.,) | in order for the number of 
overshoots of the saturated step response to remain the 
same. 

Although the above analysis has been rather crude, 
it does indicate the general effect of a decrease in the 
rate loop gain crossover frequency. If more detailed and 
accurate information is desired, it is suggested that 
analog computer simulation be employed, because an 
extension of the analysis is quite cumbersome. 


An 
Vin 


We 


We ~ 


(87) 
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APPENDIX II] 
ANALYSIS OF LEAD NETWORK SYSTEM 


The lead network system was shown previously in 
Fig. 5(a). This Appendix shows that during the satu- 
rated step response the lead network system of Fig. 
5(a) behaves in the same manner as the equivalent rate- 
feedback system of Fig. 5(b). 

Fig. 5(a) shows that the acceleration A, of the lead 
network system is related to the error X, by 


Ag = wels + w1)Xe = we(Xe + wiX.), (88) 


where X, is the error rate, that is, the time derivative of 
the error X,. Since X, is equal to X,—X-,, the error rate 
is equal to 


Veneer ec Gay od (89) 
Substituting (89) into (88) gives 
Ales = we|wiX = V.| —+- Ee am (90) 


At the instant the step is applied, X, is infinite; after 
that, X, is zero. Thus, X, is an impulse of area equal to 
the size of the step input. If this impulse is neglected, 
A, is equal to 


A, = w,|wiXe — Vel. (91) 


This same equation applies to the tachometer system in 
Fig. 5(b). Thus, if the impulse term is neglected, the 
lead network compensated system of Fig. 5(a) behaves 
in the same manner as the analogous rate feedback 
compensated system in Fig. 5(b). 

In the analogous rate feed system, the position loop 
gain crossover frequency is equal to the lead network 
break-frequency w, of the lead network system, and the 
rate loop gain crossover frequency is equal to the gain 
crossover frequency w, of the lead network system. 
Hence the previous equations that were derived for the 
rate feedback system apply to the lead network system 
if w, is replaced by w;. Thus, for the lead network sys- 
tem, the lead network break-frequency w, is related to 
the number of overshoots of saturated step response by 


2, A m. 
J i 


w) = (1+ 2). (92) 

Now consider the X, impulse that was neglected in 
(90). The lead network transfer function in Fig. 5(a) 
was approximated as (s+a;), but, as was shown by (24), 


is actually 
Ss — WO] 
(s/aw,) + 1 


(This transfer function includes an amplifier gain factor 
of a.) If the effect of the lead network factor [(s/aw,) +1] 
is considered, the term (w,X,) in (90) becomes 


Oe ( SA ) 
———— = w,aw, | ——— }. 
1+ (s/aw7) s+ aw; 


Designate the inverse transform of this term as AA,. 
For a step input of magnitude X, the transform of X; is 
equal to X/s, and the inverse transform of (55) is 


(93) 


(94) 


DN Ah, == MA 


(95) 


Assume that the step input is large enough for the 
system to accelerate to maximum velocity. Designate 
as ft, the time for the system to accelerate to maximum 
velocity [#.e., (tz—f) in Fig. 3]. The time é. is equal to 
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iV Gn A ae (96) 


During the time ¢, that the system is accelerating, the 
exponential term of (95) is of no concern because the 
motor is completely saturated anyway. This exponential 
merely tends to drive the motor more into saturation. 
At the end of this period, the term is equal to 


AA.(ta) = (Xw,aw7) eta = (Xeweaw))e~wVmlAm, (97) 
Substituting (92) into (97) gives 
AAs) = (Xo,ao)er2, (98) 


A reasonable value for a is 10. For no overshoots (7.e., 
for n=0), the exponential factor of (98) is e~®°, which is 
equal to 210°. This factor is even smaller when u is 
greater than zero. Thus, the exponential term is essen- 
tially zero at the time ¢, (when the motor has reached 
maximum speed). 

The above discussion has justified the approximation 
of neglecting the X, term in (90) when calculating the 
saturated step response of the lead network system. 
The transient due to the X, term is of no consequence 
provided that the input step is sufficiently large to ac- 
celerate the motor to maximum speed. By the time the 
motor has come out of saturation, this transient has 
vanished. 


Analog Study of Dead-Beat Posicast Control” 


G. H. TALLMAN} anv O. J. M. SMITHt+ 


Summary—A method is presented which eliminates the oscilla- 
tions and overshoot in a lightly damped servomechanism in a time 
of considerably less than one cycle of the uncompensated oscillation. 
The frequency response of a resonant system compensated in this 
manner can be made flat up through and beyond the resonant fre- 
quency of the system. For dead-beat response in one-half period, 
the frequency response is down only 3 db at the resonant frequency. 
Systems designed to have dead-beat responses in much less than 
one-half period have correspondingly wider bandwidths. 

The posicast method involves splitting up the input excitation into 
several fragments. Each fragment is applied at such a time that the 
sum of all the transient terms is equal to zero after the last excita- 
tion. The method is general and can be applied to electrical, mechani- 
cal, and even pneumatic systems. 


INTRODUCTION 


F a linear servomechanism with very lightly damped 
poles is excited with a step input, the usual result is 
many cycles of oscillation with an eventual approach 


* Manuscript received by the PGAC, July 20, 1957. 
t University of California, Berkeley, Calif. 


to steady state. If the system is excited with only half 
the original excitation, the output will display the same 
cycles of oscillation of half amplitude before reaching the 
steady-state value of one half. Suppose, however, that 
the system is allowed to oscillate for one-half cycle, at 
the end of which time itis excited with a second step of 
one half. The two transients will be 180° out of phase 
and the algebraic sum will be equal to zero. However, 
the steady-state responses to each of the inputs add di- 
rectly, and the net result is that the output reaches 
steady state in one-half cycle of the uncompensated 
oscillation. 

From the above, we can easily derive the transfer 


function of the compensator to be placed between the 


input and the lightly damped system: 
G(s) = (1+ P) =k-+ [1 — klew7/2, (1) 


k delivers a portion (or the first fragment) of the in- 
put excitation. The second portion, of magnitude (1—2), 


q 


| 
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arrives after a pure time wait of 7/2 seconds, where T 
is the period of the oscillatory system. These two por- 
tions of the excitation are spaced one-half cycle apart, 
and the transient oscillations excited by the first por- 
tion are precisely and completely cancelled by the 
transient due to the second portion if R and 1—k have 
a ratio in accordance with the damping of the transient 
oscillation. 

A comprehensive treatment of this type of control 
is given.! P is defined by (1) as (G—1). & is a constant 
that always lies between one half and one. A very 
slightly damped system will have a k slightly greater 
than one half, and a heavily damped system will have a 
Rk near unity. 


k-! = 1+ exp (aT/2) (2) 


where a@ is the damping of the poles and where T equals 
the uncompensated oscillatory period. In a stable sys- 
tem, where a is negative, 1<k“'<2. The step response 
of this compensating network is shown in Fig. 1(a) as 
an initial step to which is added another step delayed by 
T/2 seconds. If the pole-zero diagram of (1) is plotted, 
we find that the poles are all in a vertical column at 
negative infinity in the s plane (practically at very high 
frequencies) while the zeros lie equally spaced on a line 
parallel to the imaginary axis. 

If we set (1) equal to zero to solve for the positions of 
the zeros, one has the following: 


k+ (1 — k)es7? = 0, (3) 
Separating the real and imaginary parts of the exponen- 
tial with s=a+jw, 


62008. 01 f 2. = (4) 


je-*T? sin wT/2 = 0. ey 

The imaginary part equals zero when sin w7/2 equals 
zero, so that 

|cos.aF/2 | =(1.0. (6) 


In (4), e-#7/? is always positive and k/(1—) is positive 
for 0<#<1.0. Therefore, cos w7'/2 must be negative. 


cos wI'/2 = —1= (—1)* (7) 


where is odd. 
wT /2 = nr (8) 


where 7 is odd. 
The vertical positions of the zeros in the s plane a are 


therefore at 
w = n(27/T). (9) 


a 10, J. M. Smith, “Posicast control of damped oscillatory sys- 
tems,” Proc. IRE, vol. 45, pp. 1249-1255; September, 1957. 
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(a) 


(b) 


Fig. 1—Step response of compensator for complex 
and real s-plane poles. 


From (4) and (7) 
eet? = k/(1 — k) (10) 
exp (eT /2) =k +1. (11) 


From the above, (2) can be derived. All of the zeros of G 
have the same @ coordinate in the s plane as the a of 
the system poles being cancelled. 

Posicast compensation can then be described as a 
process whereby complex zeros are generated so as to 
fall upon the complex poles of the system, thereby mak- 
ing transient oscillations impossible. A single step input 
to the complete system produces only a single pole at 
the origin. Two steps, however, appear at the input to 
the resonant component. After the first step, it seems 
as though the zeros are not present, and the system will 
start an oscillation. After the second step, which occurs 
a half cycle later, the second component of oscillation 
cancels the first. 

An alternate, but equivalent way of looking at posi- 
cast control is to consider the component oscillations as 
phasors in time phase. The two oscillations are repre- 
sented by two phasors displaced 180°. 

After the first excitation, we have one phasor repre- 
senting the first component of oscillation. After the sec- 
ond excitation, we have two phasors 180° apart. These 
phasors will cancel after the second excitation if the 
magnitudes are equal at t=T7/2. To achieve this, the 
magnitude of the first excitation must be larger because 
of attenuation due to the damping. 

The phasors representing the two components of os- 
cillation are shown in Fig. 2. The first component is 
shown as decreasing in magnitude as it rotates counter 
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Fig. 2—Phasor diagram for the oscillatory components excited by 
the double-step output of the compensator. AC vectors. 
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Fig. 3—Phasor diagram for the steady-state zero-frequency com- 
ponents excited by the double-step output of the compensator. 
DC vectors. 


clockwise with increasing time. When the first phasor 
has rotated 180°, the second component is introduced to 
cancel the first. 

The dc components of the response are shown in Fig. 
3. However, the de phasors do not rotate in time and, 
therefore, the steady-state output is given by the sum of 
the two dc phasor lengths. 


APPLICATION TO SYSTEMS 


Fig. 4(a) shows an oscillatory negative feedback sys- 
tem. Fig. 4(b) shows the compensator which is placed 
ahead of the system input. Ko is k in our previous equa- 
tions. K; is (1—k). Fig. 4(c) is the step response of the 
compensator and system combination. 

The derivative of the step response shown in Fig. 
4(c) is the impulse response of the combination of com- 
pensator and resonant component. 

This system is linear, and its transient response to 
any arbitrary input is the convolution of the derivative 
of the response shown in Fig. 4 and the input function. 
For all possible inputs, the output arrives at its final 
value tangentially in a time 7/2 after the input function 
stops varying and remains constant thereafter at its 
final value. 

The transference of (1) was realized by using a trans- 
mission line, two potentiometers, and a summing juncs 
tion, as shown in Fig. 4(b). Other methods of using a 
transmission ‘line or its low-frequency analog can be 
used, but the method shown here is especially appli- 
cable to analog computers. 

The two-pole system is the simplest case for which 
compensation is necessary to remove a resonant peak 
and, therefore, was the first case treated. The system 
constructed was that shown in Fig. 4(a), consisting of 
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Fig. 4—(a) Two-pole system, (b) two-pole system compensator, 
(c) output for step input. 
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Fig. 5—Open-loop control. 


negative feedback around a cascade of a dc gain, a time 
constant, and an integrator. This gave an open-loop 
s-plane configuration of two poles and no finite zeros. 
Our hypothesis was that the use of the compensator in 
Fig. 4(b) ahead of the system in Fig. 4(a) would result in 
a step response for the entire system like that calculated 
and plotted in Fig. 4(c). 


COMPUTER RESULTS 


An artificial delay line was constructed by a cascade 
of constant R sections, having a total delay of approxi- 
mately 3.25 msec. There have been many articles 
recently on better approximations to a distributed 
parameter delay line,?~“ but since the system considered 
was a low-pass filter, small deviations from a pure dis- 
tortionless delay were not noticed in the output. The 
system which was constructed on the computer to 
match the block diagram of Fig. 4(a) is shown in Fig. 
5. The oscillation frequency could be varied from 40 
cps to approximately 1000 cps. 

The frequency was set at 154 cps with the values 
shown in Fig. 5. If the frequency of oscillation of the 


2S. Darlington, “Network synthesis using Tchebycheff poly- 
nomial series,” Bell Telephone Sys. Mono. no. 2003; July, 1952. 

§L. Storch, “Synthesis of constant time delay ladder networks 
using Bessel polynomials,” Proc. IRE, vol. 42, pp. 1666-1675; No- 
vember, 1954. 

4E. S. Kuh, “Synthesis of lumped parameter delay lines,” 1957 
IRE CoNnvVENTION RECORD, pt. 2, pp. 160-174. 
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analog simulation in Fig. 5 is computed from element 
values, it should be 92 cps. However, the capacitors 
used in the integrators were of a common paper and foil 
type, whose nominal values were not checked and prob- 
ably did not have the values shown. It was assumed 
that the individual computing amplifier gains were 
each infinite. The actual gains were not measured. The 
frequency of 154 cps was stable and easily reproducible 
over a period of several months. 

The dimensionless damping, £, was set at approxi- 
mately 0.1 to give a highly oscillatory transient. A 
smaller damping ratio would have been desirable, but 
it was found that with too small a damping ratio the 
pole locations would drift, sometimes into the right 
half s plane during the operating time. This was prob- 
ably due to slow changes in amplifier gain as a function 
of temperature. With this circuit, adjustment of either 
the time constant or the dec gain varied both the fre- 
quency and damping ratio. 

The compensator is shown in Fig. 5 to consist of the 
delay line with two inputs toa summing amplifier. The 
damping ratio of the compensator zeros is controlled by 
varying the input potentiometer. The two 100-kQ re- 
sistors in series with the potentiometer arms were added 
in order to get finer control on the adjustment of the 
damping ratio. 

The step response of the uncompensated system is 
shown in Fig. 6(a). The step response of the compen- 
sated system is shown in Fig. 6(b). The effect of the 
compensator is spectacularly apparent. Notice that the 
transient is completely cancelled in less than one cycle 
of oscillation. There is only one discrepancy in the re- 
sponse which would not be predicted by the above anal- 
ysis. There is a slight discontinuity in slope at about 
three-quarters amplitude up the leading edge of the 
output. The reason for this is shown in the step response 
for the delay line alone in Fig. 6(c). 

Notice that there is also a slight additional roughness 
in the line response alone after the delay time. This 
roughness could introduce additional transients into 
the system so that the transient sum is not equal to zero 
after the last excitation. By slightly altering the ratio of 
the line delay to the system half period, however, the 
sum of the transient components can be made exactly 
zero as shown in Fig. 6(b). 

The effect of a slight misalignment of the zeros is 
shown in Fig. 6(d) and 6(e). In Fig. 6(d) the resonant 
frequency, or w coordinate of the zeros, is reduced so 
that they move together along the straight line directly 
between the poles. The result is an increase in the 
residual oscillation, the magnitude being proportional 
to the s-plane distance between the zero and the pole. 
The wave is a negative sine wave after the second pulse. 
The phase of the error is +90°, due to the angle con- 
tribution of the zero nearest the pole. This was achieved 
by making the timing between the pulses too great, so 
that the second pulse arrived too late. 
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(b) 


Fig. 6—(a) Uncompensated system step response, (b) compensated 
system step response, (c) delay line step response, (d) effect of 
incorrect adjustment of the compensator. The zeros are at the 
correct damping but too low a frequency. The transmission line 
was too long. (e) Effect of incorrect adjustment of the compensa- 
tor. The zeros are at the correct frequency but at too great a 
damping. The transmission line delayed step was too small. The 
zero misalignment in this case is 90° from the misalignment in (d). 


In Fig. 6(e) the zeros were adjusted to the correct 
frequency but were shifted along a line of constant 
frequency in the s plane to too great a value of damping. 
The zeros were to the left of the poles. The result is 
again an error oscillation. The magnitude was propor- 
tional to the distance between the zero and the pole. 
The wave was a positive cosine wave after the second 
pulse. The phase of the error was zero degrees, which 
was the angle contribution of the zero nearest the pole. 
This was achieved by making the magnitude of the sec- 
ond pulse too small. 


FEEDBACK CONTROL 


It is evident from Fig. 4 that the method outlined 
above is valid for disturbances occurring at the input 
only and passing through the compensator. If a small 
disturbance occurs at any other point in the system 
the system will oscillate as if there were no compensa- 
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tion. To remedy this situation, the delay line can be im- 
bedded into the feedback loop so as to correct for dis- 
turbances anywhere in the forward channel. The result 
of this process yields the block diagram of Fig. 7.°° The 
pulse generator, P, which is shown in the diagram, is 
derived from (1) by subtracting a unit step at the origin. 
The transfer function of the pulse generator is then: 


(ae the ee yer Op (12) 


This generator is easily realized by the same network as 
that used in Fig. 4 for open-loop control, except that the 
generator delivers a step output of —[1—A| initially 
and then follows this with a step of height +[1—k]| 
after T7/2 seconds. The function P has no steady-state 
gain. 

In order to realize the closed-loop control, it is neces- 
sary to generate an approximation to the inverse of the 
open-loop function. Since the open-loop function con- 
tains two poles, one at the origin, the inverse will con- 
tain two zeros, one at the origin. Since physical realiza- 
bility dictates that this network will have some poles, 
however, an extra zero is introduced. This zero is made 
variable so that it can be used to cancel the effect of the 
poles at the operating frequency. The basic requirement 
of the network is that it have 180° phase shift at the 
operating frequency. The realization is shown in Fig. 8 
to consist of two zero-producing sections consisting of 
shunt RC sections followed by a series capacitance to 
provide a zero at the origin. It is desirable to place the 
poles of this network as far out from the origin as possi- 
ble in order to make the phase lag of the poles negligible 
at the operating frequency. As the poles are shifted out- 
ward, however, the gain of the circuit goes down, and so 
amplifiers are required to make up for the attenuation. 
The poles were placed as far out as possible while keep- 
ing the output signal above the noise level. The unit was 
then heavily shielded to reduce pickup. 

The actual computer wiring is shown in Fig. 9, oppo- 
site. The system was the same as that used for the pre- 
viously described tests. The circuit shows the additional 
elements required in this case added as a minor feedback 
loop consisting of the network for generating 1/G,G, fol- 
lowed by the pulse generator with positive feedback. 
This system was also adjusted for a dead-beat response. 
These adjustments were made in three steps. First, the 
system model was adjusted by inserting a vacuum tube 
voltmeter at the input to the pulse generator or the out- 
put of the adder from the feedback and the network for 
the model of the reciprocal function. The surplus zero 
position was adjusted until there was a null or mini- 
mum of transmission at 154 cps. Next, the gain of the 
branch paralleling the model of the reciprocal system 
was adjusted so the magnitude of the null corre- 
sponded to a complex zero whose Q was essentially equal 


5 See Smith, op. cit., for more details. 

§O, M. J. Smith, “Mixed distributed and lumped parameter sys- 
aoe oe IRE WESCON Convention Record, pt. 2, vol. 1, pp. 
122-132. 
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Fig. 7—Closed-loop control. 
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Fig. 8—Simulation of 1/G,Go. 


to the closed loop Q of the uncompensated system. 
Lastly, the voltmeter was removed and the system ex- 
cited with a square wave of about 20 cps. The feedback 
resistor around the output amplifier of the pulse gen- 
erator was then adjusted to minimize the transient. The 
transient was not eliminated entirely but was signifi- 
cantly reduced in amplitude. The error was probably 
the result of the combination of the error in the delay 
line and the error in the network performing the inverse 
of the open-loop function. As in the previous case, this 
error was not important because the ratio of the trans- 
mission line delay time to the system period could be 
adjusted so that all transient inputs would have a 
dead-beat response. 


More CoMPLEX SYSTEMS 


In Fig. 1 there are drawings of two different step re- 
sponses of posicast compensators. The second step re- 
sponse defines the compensator which should be used for 
a real pole. If the transcendental equation of (1) is now 
examined for poles and zeros with k>1, it is found that 
the poles have the same location as before, but the value 
of 2 in (7)—(9) is even instead of odd. The zeros have 
the same value for a but are located at frequencies of 
n(2r/T) where nm is even. There is one zero on the 
real axis of the s plane and a column of complex zeros. 
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Fig. 9—Closed-loop control. 


Notice that it takes only one degree of freedom to de- 
termine the real zero, and the other degree, the time 
delay, might be used for other poles. In other words, the 
time of cancellation for a real pole can be anything de- 
sired within the limits of physical realizability and the 
maximum corrective signal. 

The method of compensating more complex systems 
is now clear. We assemble a compensator for each pole, 
or complex pair of poles, and cascade the entire system, 
giving a transference: 


H= et (ks + [1 — k;Je*7). (13) 


i=l 
A simpler form to realize becomes evident, however, if 
we convolve each of the above transfer functions to- 
gether: 


HS Ree (14) 


t=1 

where a 7; is required for each pair of resonant poles and 
for each real pole. This function can be realized by the 
same method that has been used thus far, 7.e., by taking 
taps off a delay line at the T; and applying these points, 
multiplied by the proper constants, to a summing junc- 
tion. For two sets of complex poles, (14) simplifies down 
to four terms: 


H = By a= Boe8?1 aL B3e 87? + Bae 8 Prt 72], (15) 


The zeros in this expression seem to go through an in- 


teresting sequence in which they are printed into the 
s plane at the time 7; to compensate for the highest 


frequency poles, removed from these poles at time 7», 


ae 


and finally applied again at time 7:+7>. At this last 
time the zeros are on top of both sets of poles. Since a 
real pole compensation does not depend on the time in- 
volved, the above expression will simplify even further 
for three-pole systems. 


FASTER RESPONSE TIMES 


The above analysis is based entirely on the desira- 
bility of one-half-cycle response. The control method is 
general, however, and easily affords transient times of 
less than one-half cycle. To see the method for this type 
of response, we make the substitutions: 


g = ef 


(16) 


where J; is the time spacing between the two pulses 
closest together in the compensator and is very much 
less than the period T of the uncompensated system. 


Pier. (17) 


This transforms the jw axis of the s plane into the unit 
circle in the z plane with other lines of constant a being 
mapped into concentric circles about the origin. The 
zeros of the delay line compensator are now easily 
plotted onto the z plane. For the half-cycle response 
when 7;.= 7'/2, the z-plane configuration consists of one 
zero on the negative real axis. The position of the zeros 
in (1) is determined by the following: 


z= — [1 —2]/e. 


If 7, is set equal to 7/4, there is a change of a factor of 
two in the base for the zg plane (basing the plane on one- 
quarter cycle rather than one-half cycle) and the zeros 
of (1) will be on the imaginary axis of the z plane. 


(18) 
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2= + jV(1 — k)/k. (19) 


If transient times of less than one-half cycle are de- 
sired, it is necessary to move the zeros on a circular arc 
towards the positive z axis. By this means, since the fre- 
quency of the lowest zero is fixed by the system to be 
compensated, the next zeros in the columnin the s plane 
are moved to a much higher frequency. This results in 
much greater bandwidth, and consequently shorter 
transient times. The transient time is 27’,. Note that the 
base of the z plane must shift with the zeros so that the 
radian frequency of the lowest order s-plane zero will re- 
main the same at w of 27/7. The perunit system tran- 
sient time of 27T;,/T in radians is given by the total-in- 
cluded angle between the zeros in the z plane. This 
angle, which is shown as 20 in Fig. 10, is equal to the 
number of degrees of electrical oscillation of the entire 


system for a step or impulse input. 
0 = QnT,/T. (20) 


The transfer function of the compensator for these re- 
sponse times can be computed from the gz plane to be: 


H = [B22 + Boz + Bslz-2. (21) 


The zeros of the quadratic expression are then found 
and the corresponding coefficients are readily evaluated. 

For no damping, the zeros lie on the unit circle in the 
z plane, and the coefficients in (21) are: 


Bi = Bs = 1/[2 — 2 cos 6] 
B = — [cos 6|/[1 — cos 4]. 
6 is the half-transient time in radians of the uncompen- 


sated oscillation. An example of this method is given by 
the coefficients for one-third cycle-response: 


iPaee pave otis e386 


(22) 
(23) 


(24) 


where 7 equals the uncompensated oscillation period. 

As the zeros are moved closer to the real axis, the 
trend of the coefficients is to get larger as the axis is 
approached. Response times less than one-sixth cycle are 
not feasible unless unusual power capabilities are pres- 
ent in the system. In the limit, the step response of the 
compensator consists of a unit doublet and a unit step, 
both at the origin of time. The method of realization of 
the transfer function of the compensator for faster re- 
sponse times than a half cycle is shown in Fig. 10. This 
compensator is also an extension upon the basic com- 
pensator developed above. 


EFFECT OF MISALIGNMENT 


The above derivations specify that the pole is exactly 
cancelled by the zero. For some special purposes it is 
advisable to introduce a slight misalignment. For exam- 
ple, a slightly greater bandwidth can be obtained by 
moving the zero up and to the left of the pole. This in- 
crease in bandwidth is bought at the expense of a slight 
overshoot in the step response. 

For sinusoidal operation at the transient resonant fre- 
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Fig. 10—(a) Fractional cycle transient zeros, 
(b) complex zero generator. 


quency, the system sinusoidal gain when the zero lies 
exactly on top of the pole is approximately 7/4. The 
exact gain, when the zero has the same w value as the 
pole but has a different a value, is 


r [ 1 ete Ge) | 
4 Lata/1 — (n/2)?1 


¢, is the natural transient dimensionless damping ratio 
of the pole of a/w, or tan sin~! (a/wo). § is the damping 
ratio based on the s plane radial distance wo. € is a/wo or 
sin tan—! ¢,. From these relationships it can be shown 
that 


(25) 


Gamaiaw— 
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(26) 


6 is the ratio of the a value of the zero divided by the a 
value of the pole. It is the perunit error in the value of 
either a or ¢, for the zero, measured positively to the 
left in the s plane. 

For small values of 6 and of ¢., an approximate ex- 
pression for the system gain at the transient frequency is 


1+ 61 — n) — 42€n 
( Tn) a], (27) 
1 — ¢ 


Tv 
Gain at wo, = =| 
4 
or very approximately, 
: Tv 
Gain & ve [1 + gl. (28) 


A good approximation for the maximum sinusoidal 
resonance for a wide range of values of 8 is 


1 


(+. + a a1 — ¢) 


Moe ar (29) 
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Either (25) or (29) can be used to solve for the desired 
value of 8, given the maximum resonant rise, Minax. 
This gives the amount of zero-location shift from the 
best transient response (dead-beat) as a function of 
bandwidth specifications. 

In addition to desirable misalignments, there can be 
undesirable drifts in the relative positions of the zeros 
and the poles. Closed-loop feedback systems that con- 
tain machinery or nonlinear devices may have poles 
whose positions move as a function of signal level. Am- 
plifiers and other electronic components may change 
gain and transfer functions with temperature changes. 
The net result of these effects can be represented as a 
pole motion, and (25), (28), or (29) give the maximum 
change in JM at the sinusoidal resonant frequency as a 
function of the perunit damping misalignment 6 of the 
zero with respect to the pole. 
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CONCLUSION 


A method of control has been developed which brings 
transients to their steady-state values tangentially in a 
small fraction of the normal oscillatory period by cancel- 
ling the s-plane resonant poles with z-plane zeros. The 
method gives excellent waveform. The frequency re- 
sponse of the compensated system is down only 3 db at 
a frequency equal to 47, where J is the arbitrarily 
chosen transient time. The method is easy to apply, the 
only difficult hardware being an artificial transmission 
line. 

This mode of operation can be applied to a system 
with an arbitrary number of poles by a simple extension 
of the two-pole case which has been discussed here. The 
steady-state output can be obtained in as short a time as 
desired, the only restrictions being the dynamic range 
and snr of the system. 


On Closed-Form Expressions for Mean Squares 


in Discrete-Continuous Systems” 
JACK SKLANSKYt 


Summary—When a system is to be optimized with respect to 
the mean square of some variable, a closed-form expression for 
that mean square is usually desired. The problem of obtaining such 
expressions for discrete-continuous systems—i.e., systems made up 
of both sampled-data and continuous subsystems—has been a 
difficulty in the past. The reason for this is that the spectral densities 
of the variables of interest often contain rational functions of 
exp (j27fT) combined multiplicatively with rational functions of f, f 
being the frequency coordinate of the spectral densities, and T the 
sampling period. Presented here is a technique for finding the de- 
sired closed-form expressions. It is based on the relation 


[ Prenowas = ¢ P*(e)Q*(@)eds, 


joo 


where Q* (z) is the ‘‘Z-transform” of Q (s), 

To illustrate the technique, closed-form formulas for the output 
and ripple of discrete-continuous systems and for the control error 
of sampled-data feedback systems are derived, and an application to 
a ‘“track-while-scan” system is given. 


I. INTRODUCTION 
IVE con: “sampled-data” systems in the fields of 


communication and control are “discrete-con- 
tinuous”; .e., they are combinations of both 
purely “discrete” and purely “continuous” subsystems, 
the terms “discrete” and “continuous” referring to time- 


* Manuscript received by the PGAC, December 6, 1957. 
+ RCA Labs., Princeton, N. J. 


domain characteristics of the signals involved. To evalu- 
ate the performance of such systems, it is often desired 
to compute the mean square “system error”—i.e., the 
average of the squared difference between the desired 
output and the actual output—in response to a random 
input. 

Furthermore, it often is desired to optimize the sys- 
tem with respect to one or more dynamic parameters. 
To do this conveniently, a closed-form expression for the 
mean square error in terms of the dynamic parameters 
of the system is required. In addition, closed-form ex- 
pressions for other signals in the system such as the 
output, or related functions such as the ripple may be 
needed. } 

In the past, these mean squares have been expressed 
in integral forms requiring numerical methods for their 
evaluation.!~4 When closed-form expressions were re- 
quired, they were obtained by approximating the dis- 
crete-continuous system by an all-discrete or an all- 

1 J. Sklansky and J. R. Ragazzini, “Analysis of errors in sampled- 
data feedback systems,” Trans. AIEE, vol. 74, pt. II, pp. 65-71; 
May, 1955. 

2 J. R. Ragazzini and L. A. Zadeh, “Analysis of sampled-data sys- 
tems,” Trans. AIEE, vol. 71, pt. Il, pp. 225-234; November, 1952. 

3M. Mori, “Statistical treatment of sampled-data control systems 
for actual random inputs,” Trans. ASME, vol. 80, pp. 444-456; 
February, 1958. ; ; 

4R. M. Stewart, “Statistical design and evaluation of filters for 


the restoration of sampled data,” Proc. IRE, vol. 44, pp. 253-257; 
February, 1956. 
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continuous system. Thus, a need for closed-form, ana- 
lytical evaluations without approximations has been 
apparent. This paper answers this need for a class of 
situations frequently encountered in practice. 

Given in Section II is a formula for the mean square 
output for systems in which: 1) the system transfer 
function can be expressed as a product-of a rational, 
physically realizable transfer function of the Laplace 
operator s and a rational, physically realizable transfer 
function of e’? (or a sum of such products), 2) the in- 
put’s spectral density is a known rational function of 
frequency, and 3) the input signal is independent of the 
sampling process. In addition, of course, the over-all sys- 
tem must be stable, since otherwise the mean squares 
would be infinite, and the computed values of the mean 
squares, if finite, would be incorrect. 

The techniques used in deriving this formula can also 
be applied to finding closed-form expressions for other 
time functions related to the system, such as the system 
error, the control error,5 and the ripple. These tech- 
niques are described in Appendix II. Formulas for the 
control error and the ripple for error-sampled feedback 
control systems are developed in Appendix ITI, making 
use of the techniques of Appendix II. 

At this point an equation displaying the essential 
mathematical contribution of this paper is presented. 
This contribution amounts to a method of evaluating 
the integral over all real frequencies of a product of two 
spectral densities, one of which is rational in e*” and the 
other rational in s. This integral is the left member of 
the equation, given below, and the method of evaluating 
it is indicated by the right member. The basis for this 
relationship is derived in Appendix II. 


fs P* (eT) -Q(s)ds 


—jo 


2Qrj 
1 
ee —$ P*(z)Q*(2)z-'dz, (1) 
TT) unit circle 
where 


O* (er?) = Z0(s) 
O(p) 


Qrj c—jo 1 Cue (s—p) TP 


Q*(z) 
K 1 c+jo 


dp. (1a) 


It is assumed in (1a) that ¢ is chosen large enough to 
cause the line of integration, c—j~ toc+j~, to be to 
the right of all the poles of Q()). 


II. AN ANALYTICALLY-EVALUABLE FORMULA FOR 
THE MEAN SQUARE OUTPUT 


The discrete-continuous system illustrated in Fig. 1 


5 A distinction is made here between the system error and the 
control error. The former is the difference between the desired output 
and the actual output; the latter is the difference between the input 
of a feedback system and the feedback signal. 

6 “A” means “equals by definition.” Derivation of (1a) given by 
E. I. Jury, “Analysis and synthesis of sampled-data control sys- 
tems,” Trans. ATEE, vol. 73, pt. I, pp. 332-346; September, 1954. 
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is considered. The input is a stationary random signal, 
u(t). This is sampled at the rate 1/T to produce the 
train of impulses 


io] 


ut(t) = SS u(nT)a(t — nT). (2) 
A linear system, whose transfer function is the product 
of A*(es?) and H(s), responds to this “sampled” input 
to produce the output v(#). It is assumed that A*(e*7) 
and H(s) are both physically realizable transfer func- 
tions, that the former is rational in e*7, and that the 
latter is rational in s.7 In addition, it is assumed that 
v(t) remains bounded for all ¢. This latter assumption 
will hold good when A*(z) and H(s) are “stable,” 7.e., 
when the poles of A*(z) lie inside the unit circle of the 
z plane and the poles of H(s) lie inside the left half of 
the s plane.® 


K(s) 
(ae eee 
ie) 7 eet x ost v(t) 
uth) A (e®') H(s) 
INPUT T OUTPUT 


Fig. 1—A common type of discrete-continuous system. 


Under these assumptions, w(t) and v(t) are related by 
V(s) = U*(e7?) A*(e°7) H(s), (3) 


where U*(e’”) is the z transform of u(t). (“gz transform” 
is just a special designation for the Laplace transform 
of (2), with zg replacing the variable e*?.) Thus 


U*(s) = YO ulnT)s. (4) 
Techniques of handling transforms of this type are dis- 
cussed.!: The latter reference is of particular interest 
here, because it includes some discussion of bilateral 
z transforms, which are associated with the autocorrela- 
tion functions of sampled random signals. 

Assuming that the input’s autocorrelation function, 
Y.(r), is known, and that the corresponding spectral 
density, V.(jw), is a rational function of w, the following 
equation is an analytically evaluable formula for the 
mean square output of the system in Fig. 1. [For con- 
venience we use the symbolic form Y(jw) to denote spec- 
tral density rather than the more usual notations V(w) 
or V(w?).] 


weve i d 
A = ft Wu¥(z) A*(2) AM(s We) (5) 
2rjJ unit circle Tz 


* Convenient, necessary, and sufficient conditions for this assump- 
tion are given in A ppencee Le 

8 J. G. Truxal, “Automatic Feedback Control System Synthesis,” 
McGraw-Hill Book Co., Inc., New York, N. Y.; 1955. 

* W. K. Linvill, “System theory as an extension of circuit theory,” 
He TRANS. ON Circuit THEORY, vol. 3, pp. 217-223; December, 

0 The fraction dz/(72) is written as shown, because dz/(Tz) =ds, 
thus exhibiting the equivalence of (5) to the integral of a spectral 
density over the real frequency axis. 


lee tea ls 


+ 


— - 
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where 
Mee) = Sy ValnT) se, (6) 
Wi*(z) = DS valnT)2-. (7) 


n=—oo 


The function ¥,(zT) is defined as the autocorrelation 
function of £~! H(s), with the correlation interval equal 
to nT. Thus, 


Vilt) = if h(i)h(t + z)dt, (8) 
where h(t) is the inverse Laplace transform of H(s), i.e., 
h(t) is the impulsive response of the continuous portion 
of the system of Fig. 1. 

The derivation of (5) is given in Appendix II. Tech- 
niques of evaluating it analytically are discussed in 
Section ITI. 

Although (5) explicitly pertains to the output, the 
techniques involved in its derivation can be used for 
deriving closed-form expressions for other mean squares 
related to the system, such as the system error, control 
error, and ripple. To illustrate this fact, analytically 
evaluable formulas for the mean square control error 
and the mean square ripple in error-sampled feedback 
systems are derived in Appendix III. 


III. EVALUATING THE FORMULA 


To evaluate (5), the closed-form expressions for 
W,,*(z) and W;,*(z) are found. Then these expressions are 
combined with A*(z) and A*(z~!) to obtain the inte- 
grand of (5), and the residues of this integrand inside 
the unit circle are obtained. 

To find W,*(z), use either (6) in conjunction with a 
table such as Barker’s!"! or use the relation 


W.*(z) = ZW.(s), (9) 


where 


Vu(s) = Wal jor) |urmet (10) 


and where Z represents the z transformation of the time 
function upon whose Fourier or Laplace transform Z 
operates, z.e., Z represents a mapping from the s to 
the z domains. A good table of such mappings is given 
by Barker;!! shorter tables are given by Truxal? and 


Jury.® An integral representation of the Z transforma- 


tion is given by (1a). 

To illustrate the method of using the tables to 
find Y*(z), consider the case where the spectral density 
of u(t) is a/(w*+a?). Replacing w? by —s’, W.(s) 


11 R, H. Barker, “The pulse transfer function and its application 
to sampling servo systems,” Proc. IEE, pt. 1V, monograph no. 43; July 
15, 1952. Note: to use the tables, one would have to express Wul(r) asa 
sum of two “one-sided” time functions, one of which is zero for nega- 


~ tive time, and the other zero for positive time. The reason for this is 


that the present tables of z transforms were computed only for one- 


, sided time functions. 
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=a/(—s’+a?). The Z transform of this expression may 
be found in a table or from (1a), where 


104 
Op) = a 
(p — a)(p + a) 
The obtained transform is 
—z sinh aT 


W,*(z) = ; 
@) 2? — 2z coshaT + 1 


(Even though the tables list only the transforms of one- 
sided time functions, the tables can be used for two- 
sided time functions when going from the s to z or from 
the z to s domains. Essentially, this is because both the 
negative and positive sides of any complex exponential 
of the form Ki e@+8)'—which is the type of time-func- 
tion component contributing partial fractions to the 
s and z transforms—will have the same s transform, 
and both sides will also have the same z transform.) 

W,*(z) may be found in a similar manner by using (7) 
or the relation: 


W,*() = Z| H(9) H(—s)). (11) 


This relation is a consequence of (7) and the fact, im- 
plied by (8), is that the Fourier transform of y,;(r) is 
H(s) H(—s), with s replacing jw. 

To illustrate, suppose H(s) =1/s. Then, 


(12) 


IV. AN ALTERNATE EXPRESSION 


The following alternate formula for v(t) stems di- 
rectly from (5) through applying in reverse order the 
techniques used in progressing from (32) to (34) in Ap- 
pendix II. This formula is of some academic interest, 
because the integration is made over the imaginary axis 
in the s domain, as in Laplace inversion; however, it 
offers less computational convenience than (5) because 
of the inverse Z transformation required. 


1 clea | 
2(f) = — — Wx(s)ds, 13 
v°(?) rj J ye T K(s)ds (13) 
where Vx(s) is defined by 
ZYi(s) = Wy*(z) A*(z) A*(e-YWp*(s). (14) 


To find W;(s), an inverse Z transformation of (14) is 
made by expanding the right member into partial frac- 
tions with the form a;/(1—e~*i?z“!), and assigning corre- 
sponding fractions of the form a;/(s+6,;) to W;(s). 
(Here 6; and a; may be real, complex, or imaginary.) 

A computational advantage of this formula is that 
available tables! have direct applicability to its evalu- 

12H. M. James, N. B. Nichols, and R. S. Philips, “Theory of 
Servomechanisms,” M.I.T. Radiation Lab. Series, McGraw-Hill 
Book Co., Inc., New York, N. Y., vol. 25; 1947. (See table in ap- 
endix.) 
Pons R. C. Booton, Jr., M. V. Mathews, and W. W. Seifert, “Non- 
linear Servomechanisms with Random Inputs,” M.I.T. Dynamic 


Analysis and Control Laboratory, Cambridge, Mass., Rept. No. 70; 
August 20, 1953. (See table in appendix.) 
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ation. However, the same can be said of the y-domain 
version of (5) obtained by the transformation z= (1+ 4) 
/(1—y), and, except in the simplest cases, this trans- 
formation seems to involve less labor than the Z inver- 
sion of (14). Furthermore, tables can be constructed 
which would apply directly to (5) without a z-to-y trans- 
formation. 


V. ILLUSTRATIVE EXAMPLE 


The use of (5) will be illustrated by the problem of 
optimizing a particular track-while-scan system." 

Consider the sampled-data feedback system shown in 
Fig. 2. This forms a subsystem handling one of the 
cartesian coordinates in a common type of track-while- 
scan system. It is desired to choose values of a and 8 
which in some sense will provide optimum performance. 


Fig. 2—The block diagram of a common type of 
track-while-scan system. 


One of the figures of merit that has been considered for 
such an optimization is the ratio of the mean square of 
the output, v(t), to the mean square of the input, w(t), 
when u(t) is “white” noise whose frequency spectrum is 
flat over a very wide but finite band. 


IRE TRANSACTIONS ON AUTOMATIC CONTROL 


March 


(5), since the bracketed factor of (15) corresponds to 
A*(z), and the factor to the right of the brackets corre- 
sponds to H(s). 

The task now at hand is to find closed-form expres- 
sions for V,*(z) and W,*(z), which are needed for inser- 
tion in (5). To do this we shall use (6) and (11). 

Since the spectrum of u(t) is spread uniformly over a 
very large though finite band, the autocorrelation func- 
tion of u(t) is 


a for 7 = 0 (16) 
OM Seo ey eee 
Hence, by (6), 
W,*(z) = u(t). (17) 


The function H(s) H(—s) in (11) is equal to G(s) 
G(—s), where G(s) is defined in Fig. 2. Hence, by (11), 


B 
“54 

(g2 — ae 

we pay 


a? 
W),*(z) — z| - stale 
2 


(462 + 1202)z 
B? 4 6a2 


a |. (18) 


In going from the second to the third member of (18), a 
table such as Barker’s!! is convenient. (An interesting 
alternative technique for this purpose, using “modified 
z transforms,” is described by Mori.’) 

Substituting (17) and (18) into (5) and setting A*(z) 
equal to the bracketed factor of (15) obtains 


B® — 60% + (48? + 12a)z + (6? — 6a%)z? 


v*(Z) 


It can be shown that u(f) and v(t) are related by 


Vos) = va(9)| aoe G(s) 


or after making substitutions: 


ee ural 2—(2 aE coe) oe 


where z=e*?, and where the time scale has been normal- 
ized so as to have 7 =1. [A note to the reader under- 
taking to derive (15): In computing G*(z), the discon- 
tinuity in £—!G(s) should, for physical reasons, be re- 
placed by a line segment whose slope is finite, although 
large and positive. This is equivalent to subtracting the 
initial value of the time function from G*(z).| 
Comparing (15) with (3), we see that this system 


satisfies the requirements for the evaluation of v(t) by - 


4 A track-while-scan system can be broadly defined as a system 
for estimating the “present” value of a time function from a sequence 
of “past” samples of that function. 

% J. Sklansky, “Optimizing the dynamic parameters of a track- 
while-scan system,” RCA Rev., vol. 18, pp. 163-185; June, 1957. 


ADs 
aj J unit cirele O[(1 — a) — (2-—a—B)z+ 2] [1 — (2-—a—B)z+ (1 — a)2?] 


dz. (19) 


This can be evaluated either by finding the residues 
directly, or, more conveniently, by substituting 
z=(1+y)/(1—y) and using one of theavailable tables. 12:8 
The result is 


0) 
wl) 


6a? — 3a6 + 68 — B? 
3a(4 — 2a — B) 


(20) 


This is a closed-form expression that can be applied to 
choosing optimum values of a and £. (A graphical opti- 
mization scheme making use of this expression is pre- 
sented.) 


VI. CONCLUSION 


The techniques used in the derivation of (5), which 
are summarized to an important extent by (1), con- 
stitute tools for finding closed-form expressions for 
mean squares that previously could only be found nu- 
merically. These techniques make possible relatively 
straightforward optimizations of discrete-continuous 
systems with respect to selected mean squares. For 
illustrative purposes, two other formulas, namely (41) 
and (47), are derived in Appendix ITI. 
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The assumptions making these techniques practicable 
are that WV, (jw) is rational in w; that the discrete and con- 
tinuous portions of the systems are rational, physically 
realizable transfer functions in e*? and s, respectively ; 
and that the input is independent of the sampling proc- 
ess. In addition it isassumed that the over-all system is 
stable, so that the mean squares will not be infinite. 

As an aid in evaluating the formulas obtained by 
these techniques, a table of integrals around the unit 
circle for rational integrands would be convenient. Such 
a table would be similar in organization to the presently 
available tables of integrals of rational functions over 
the real and imaginary axes.!2:# 


APPENDIX [| 


PHYSICAL REALIZABILITY OF A*(e*7) AND H(s) 


To test for the physical realizability of A*(e*?) and 
H(s) when these functions are rational in e’? and_s, re- 
spectively, it is convenient to use the following necessary 
and sufficient conditions. 

Let z=e*?. Assuming A*(z) and H(s) are rational in 
z and s, respectively, they are physically realizable 
transfer functions if and only if they can be expressed in 
the forms 


TS ee Si a ig ig al 
1 = Fi biz! op bos a wes + Ona ae 


A*(z) = (21) 


and 


SS SSE ea? 


do + dis + +++ +dpys) + 5? 


where the a,’s, 0;’s, c;’s, and d,’s are all real. 

It will be noted that the physical requirement of finite 
transmission at infinite frequency |z.e., | A*(o)| <0 
and | H( 00 )| < «| is imposed by the unity term in the 
denominator of A*(z) and by the equality in the degrees 
of the numerator and of the denominator of H(s). 
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DERIVATION OF (5) 


Formulas for the spectral density of v(t) is derived. 
This part of the derivation is a modified version of 
Stewart’s.* 

The spectral density of v(¢) is, by definition, 


v,(jo) 2 lim — | Vo) P (22) 
: Boo B 

where Vz(jw) is the Fourier transform of vg(t), which is 

a section of v(t) covering a time interval B units long. 

_ The subscript B will have the same connotation when 

used with other variables. 
Referring to (3), we make the following approxima- 

tion, frequently used in derivations of this kind: 


Valju) & Us*(ei*?) A*(e*7) H(je). (23) 
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The error involved in this approximation will disappear 
presently when B is made to approach infinity. Multi- 
plying each side of (23) by its own conjugate results in 
an expression for | Va(jw) 2. When this is substituted in 
(22) we obtain 


V,(je) 


2 


1 
=( lim — | Up*(e?) ‘| A* (ei?) |2| (jw) |2. (24) 
Boo B 


Now, as shown in the literature,! Uz*(e##") can be ex- 
panded into a “Poisson sum”: 


jars 
Up*(e#?) = = Dd. Us(jw + jnwo), (25) 


nN=—o 


where wo is the angular sampling frequency. Thus the 
bracketed factor in (23) can be expressed as a sum of 
spectral densities and cross spectral densities: 


1 
lim 3 | Up*(e%7) |? 


Bow 
1 eae So amo eae 
=— ))| lim — Ua(jo + jneo)Us(jo + jmeo) (26) 
1 m,n Boo B 
where “~~” indicates the conjugate function. 


A typical term in this summation is the Fourier 
transform of 


ss Ut + r)u(t)er ment, (27) 
T? 

where the bar indicates an averaging with respect to 
time ft. Since u(t) and the sampling process are assumed 
to be independent, (27) is equal to the product of the 
time averages of u(t+7r)u(t) and e?%—™+ot, The latter 
average is zero except when =m, in which case it is 
unity. Hence every cross spectral density in the summa- 
tion of (26) is zero. The remaining terms are the “side- 
band spectra” of u*(t): 


1 
lim — | Us* (er?) |? 


Bow 


(28) 


Applying the Poisson-sum identity (25) to the right 
member of (28), we have 


1 1 

lim — | Up*(ei?) [2 = — W,*(cio? 

lim = | p¥(ei?) | 7 ete?) 

(29) 


= % [ZW,,(s) Jrncio?. 


where W,(s) is obtained from V,(jw) by replacing w? 
by —s*%. 
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Putting (29) into (24), we obtain an expression for 
the spectral density of v(¢): 


9 


il 
W, (jo) = 9 CACC? H(jw) |?. (30) 


The mean square of v(t) is just the integral of (30) 
over all real frequencies: 


9 


co ore ew eee | 
v(t) = — ( po Wet(eier) | AMC) H (jw) |?dw. (31) 


27 J _« 

This formula has presented an obstacle in previous 
investigations, because no non-numerical technique of 
evaluating it is apparent. 

To make possible a non-numerical evaluation of (31) 
we express it as an infinite sum of integrals over succes- 
sive frequency-domain periods of the starred functions: 


wo/2+nwo 1 
\f — W,*(eiT) | A*(et#7) |? | H(ju) tao} . (32) 
—wo/2+nwo ik 

Making the change of variable w+-na0—w for the nth 
interval of integration, and remembering that the 
starred functions are even functions of frequency and 
are periodic with a period of wo, we have 


v2(t) —f [ve A* (eT) [2 


2m J —w9/2 


‘eye 
2m dX | AGo + jnwo) | dw. (33) 
Using the Poisson summation rule, (25), and making the 
change of variable e’*7—z, we obtain the desired result: 


mee a 1 dz 

Ok aaa g U8 (2) A*(2) AMY 4(2) — > (34) 
20] unit circle Tz 

where 


U,*(2) = Z[H(s)H(—s)]. (35) 


Eq. (7) provides an alternative to (35) for finding 
W,*(z). The fact that (35) and (7) are consistent can be 
shown by Z-transforming the Fourier transform of (8). 

Note that the relation between the right-hand mem- 
bers of (31) and (34) is the same as between the two 
members of (1). Eq. (1) has thereby been proved, and 
may be used hereafter in the derivations of other mean 
squares. Two examples of this use are given in Ap- 
pendix ITI. 


APPENDIX III 


MEAN SQUARES OF RIPPLE AND CONTROL 
ERROR IN ERROR-SAMPLED 
FEEDBACK SYSTEMS 


With the aid of (1), or, equivalently, the technique 
used in progressing from (31) to (34) in Appendix II, 
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we shall obtain formulas yielding closed-form expres- 
sions for the mean square ripple and the mean square 
control error in feedback systems of the “error-sampled” 
type. The block diagram of such a system is shown in 
Fig. 3. 

Heuristically defined, “ripple” is the difference be- 
tween the output of a discrete-continuous system and 
some specified smoothed version of the output, as indi- 
cated in Fig. 4. (A specific definition is given.!*) “Con- 
trol error” is defined here as the difference between the 
input and the fedback signal, so that in a unity-feedback 
system it is the difference between the input and the 
output (see Fig. 4). 


Fig. 3—An error-sampled feedback system. 


DESIRED OUTPUT 


SYSTEM 
CONTROL ERROR 
ERROR 
RIPPLE 
OUTPUT 
SMOOTHED 
OUTPUT 


10) aT av i 4T 
TIME 


Fig. 4—Time-domain definitions of “ripple,” “control error,” and 
“system error” in a sampled-data, unity-feedback system. 


Using the results of Sklansky and Ragazzini,! we have 
the following formula for the mean square ripple: 


pe fo mien 
e(t) = en ie T| {eee G*(ei#T) |2 w 
_ +f? SaGe)| Ge) 
ele T? |1 +4 G(e?) dw, (36) 
where 
U,*(e7) = Z[G(s)G(—s) nee. (37) 


Again the first of these integrals can be shown to be 
equal to the mean square output, v*(t). Hence, by (5), 
the first integral of (36) is equal to 


16 Sklansky and Ragazzini, op. cit., in the section “RMS System 
Error,” use the term “system error” as an equivalent of “ripple,” be- 
cause the desired output is assumed equal to the smoothed output. 
Awareness of this will prevent confusion. 


es WOO) a 
21}. unit circle [1 + G*(z)|[1 a G*(z-})| Tz 


. (38) 


Alternatively, this expression could have been derived 
by applying (1) to the first integral of (36) after replac- 
ing eT by z. 

Again applying (1), we get the following equivalent 
expression for the second integral of (36): 


1 
ra WW ,*(z) 


pe circle [1 ts G*(z)][1 ch Gta) . Tz ; 


(39) 


TW, *(2) = Z[Vu(s)G(G(—3)]). (40) 


Replacing the integrals in (36) by (38) and (39), we get 


1 
W,*(2)W,*(2z) — — VW, *(z 
= : (z) ¥y*(z) A a” (2) 
eb) =a 
20] 


dz 
Tz 


[1+ G*(2)][1 + G@)] 


unit circle 


(41) 


This integral can be evaluated in closed-form whenever 
G(s) and W,(s) are rational functions of s. 

In a similar fashion one can find an integral yielding 
a closed-form expression for the mean square of the con- 
trol error, e(t), of the system in Fig. 3. One begins with 
a result derived by Mori:? 


1 ee 2V,(4w 


TY —x 


G(jw) 
; E Ee eat 
V,*(ei#T)W,,(jw) \ (42) 
T|1 + G*(ei?) |? 


The first term of the integrand contributes just the 
mean square input, ?(¢), to the integral. The last term 
contributes the mean square output, v?(¢), a closed-form 
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evaluation for which is provided by (38). To evaluate 
the contribution of the middle term, we write its in- 
tegral in the equivalent form, 


1 ia eee 
ere. pe ore es 


a] 
: da 
1 G* (E72) 

Since V,,(jw) =Vu(—jw), the integral of the second term 


in the brackets equals that of the first. Hence (43) 
equals 


(43) 
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Again we apply (1), this time to (44). [Reviewing the 


significance of (1): applying (1) to (44) is equivalent to 
resolving (44) into a sum of integrals over the intervals 
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and making the change of variable w+nw)—w, followed 
by e#7—z.| As a result, (44) becomes 
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where 
V,G*(2) = Z[W.(s)G(s)]. (46) 


We now replace the integrals of the first, second, and 
third term of the integrand of (42) by w2(t), (45), and 
the right member of (38), respectively. This gives us 
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which will yield a closed-form expression whenever 
W,,(s) and G(s) are rational functions of s. 
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Bibliography of Sampled-Data Control Systems 


and Z-Transform Applications’ 
H. FREEMAN} Anp 0. LOWENSCHUSSt 


Summary—tThe bibliography given here represents an attempt 
to simplify the task of searching the literature to obtain either a 
general acquaintance with the field of sampled-data control or to ob- 
tain information on certain specific aspects of it. Due to the close 
association of the Z transform and sampled-data systems, a number 
of papers dealing exclusively with the Z transform are included in 
the bibliography. 

The bibliography includes only that material which in the authors’ 
opinion represents either a significant contribution to the field or has 
tutorial value. The bibliography is arranged alphabetically according 
to the name of the author, or the first-named author in cases of co- 
authorship. A subject index with cross-references to the author list 
is also provided. 


INTRODUCTION 


AMPLED-DATA systems are systems in which 
S the data are transmitted as periodically recurring 
samples rather than as continuous data. The sam- 
pling action causes the output of such systems to behave 
in a different manner from that of otherwise equivalent 
continuous systems. Special techniques are required for 
the design of these systems and for the prediction of 
their performance. A sampled-data system can be linear 
or nonlinear, its parameters can be constant or time- 
varying, and the sample rate can be constant or variable. 
Sampled-data problems originated with the develop- 
ment of automatic-tracking radar systems. The pulsed 
nature of radar data affects the response of the tracking 
servo systems. The increasing use of digital computers 
in control system applications has awakened further in- 
terest in the sampled-data field. This is due to the fact 
that digital computers can accept inputs and deliver 
outputs only at discrete instants of time. 

The Z transform appears to provide the most direct 
method for the analysis and synthesis of sampled-data 
systems. Nearly all the recent papers in the field of 
sampled-data systems have utilized the Z transform, 
either directly or in some modified form. In addition, 
the Z transform has been applied to the numerical solu- 
tion of continuous systems, signal flow graphs, and sta- 
tistical processes. 

Relatively little of the theory of sampled-data sys- 
tems has thus far been published in book form. How- 
ever, a very extensive number of papers are available. 
The bibliography given here represents an attempt to 
simplify the task of searching the literature to obtain 
either a general acquaintance with the field or to obtain 
information on certain specific aspects of it. Due to the 
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close association of the Z transform and sampled-data 
systems, a number of papers dealing exclusively with 
the Z transform are included in the bibliography. 

Sampled-data control systems were first studied dur- 
ing World War II [36, 58], in conjunction with radar 
tracking. The discontinuous nature of radar position 
data forced the designers of control systems to surrender 
the notions of continuous systems and to examine the 
effect of sampling on control system performance. Early 
analyses were carried out largely by means of the 
Laplace and Fourier transforms. 

The Z-transform method was first introduced by 
mathematicians for solving difference equations. This 
method was applied to sampled-data control system 
analysis by Ragazzini and Zadeh [68] and Barker [4]. 
The work of Barker includes an analysis of transporta- 
tion lag (pure delay) and a comprehensive table of 
transform pairs. Both Barker and Ragazzini and Zadeh 
discuss system analysis, prediction of response, stability 
by Z-transform methods, and hold circuits. 

Digital compensators and analog networks can be 
used to improve the response of a sampled-data system. 
Barker [4] and Bergen and Ragazzini [9] discuss meth- 
ods for digital compensator synthesis, using as criterion 
the desired system response (at sample instants) to a 
given input function, which is usually a unit step or 
ramp. Franklin [30] and Bertram [10] examine this 
technique further and utilize the mean-square error 
criterion for optimizing the digital compensator. An 
alternate approach by Jury and Schroeder [41] and 
Jury [44-45] uses the modified Z transform to examine 
behavior between sampling instants and to design the 
digital compensator so as to suppress “hidden” oscilla- 
tions. 

The mean-square error criterion is also discussed by 
Sklansky and Ragazzini [76] and by Blum [11]. 

Sampled-data systems can be analyzed by means 
other than Z transforms. A comparison with amplitude 
modulation was made by Linvill [57]. Bilinear trans- 
formations have been used by Johnson, Lindorff, and 
Nordling [38, 39]. An approximation method is due to 
Brown and Murphy [19]. 

Nonlinear and time-varying sampled-data systems 
were analyzed by Kukel [50], Kalman [47], and Fried- 
land [32, 33]. Multiple-rate sampled servos are dis- 
cussed by Kranc [49]. Multipole sampled servos were 
studied by Freeman [31]. The effect of finite sample 
width has been examined by Farmanfarma [26, 27]. 
Contactor servos are discussed by Chow [25]. 
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Proposed IRE Standard Terminology for 
Feedback Control Systems 


HE following definitions have been approved only 

by subcommittees. They are extensions of the pre- 

viously approved IRE Standard terminology, but 
they are not yet IRE Standards. They are published to 
obtain comment and suggestions for changes before they 
are approved by the feedback control committee 26.0 
and the IRE Standards Committee. 

Should these definitions be approved as Standards? 
Would they be accepted as common references? Opin- 
ions will be welcomed by the editor. It is much easier to 
change proposed Standards than it is to change ap- 
proved Standards. Comments should be made as soon 
as possible. 


Acceleration Constant. The reciprocal of the accelera- 
tion error coefficient. 


Acceleration Error Coefficient. The coefficient of the 
quadratic term in the power series expansion of the 
difference transfer ratio. 


Bode Diagram. The sum of the asymptotes of the 
logarithms of the magnitude of the factors that con- 
stitute a transfer ratio plotted as a function of frequency. 


Difference Transfer Ratio. The transfer ratio of a loop 
difference signal to the corresponding loop input signal. 


Error Coefficient. A coefficient of the power series ex- 
pansion of the difference transfer ratio, when expanded 
in terms of increasing positive powers of the Laplace 
operator. 

Frequency Response. In the sinusoidal steady state, the 
complex ratio of the output phasor to its input phasor 
as a function of frequency. 

Note: If the concept of frequency response is ex- 
tended to nonlinear systems, the sinusoidal 
components being considered should be de- 
fined, as well as the conditions of test. 


Gain-Crossover Frequency. Any frequency at which the 
logarithm of the loop gain changes sign. 
Gain, Loop. The magnitude of the loop transfer ratio. 


Gain Margin. The reciprocal of loop gain at a phase 
crossover frequency. 


Note: Gain Margin is often expressed in decibels. 
Log Magnitude-Phase Diagram. A plot in rectangular 


coordinates of the logarithm of magnitude vs phase 
angle of a transfer ratio, with frequency as a parameter. 


- Magnification. A frequency response, expressed as the 
_ ratio of the magnitude of the response at any frequency 


to the magnitude of a reference frequency. 


Note: In low-pass systems, the reference frequency is 
commonly taken as zero. 


Magnification, Loop. The magnitude of the return 
transfer ratio. 


Note: This definition coincides with that of magnifica- 
tion for unit feedback, low-pass, single loop, 
system. 

This peak value of the loop magnification is a 
characteristic of the loop, indicative of loop 
stability. 


Note: 


Nyquist Diagram. A plot of the loop transfer ratio in its 
complex plane, corresponding to all real values of fre- 
quency, usually made for the purpose of examining the 
stability conditions of the loop. 


Nyquist Diagram Generalized. A plot of loop transfer 
ratio in its complex plane corresponding to some speci- 


- fied contour in the complex frequency plane. 


Operationally Stable. A system is operationally stable 
if the utilized output is under control of its input. 


Overshoot, Relative. In the response of an automatic 
feedback control system to a step function input, the 
ratio of the extreme value of the output variable minus 
the final value, to the final value minus the initial value 
of the output variable, provided the final value is con- 
stant and not infinite. 


Phase-Crossover Frequency. Any frequency at which 
the phase angle of the loop transfer ratio passes through 
180°. 

Phase Margin. In a feedback control loop, at any fre- 
quency, the angle equal to 180° minus the absolute 
value of the phase of the loop transfer ratio. 


Phase Margin at Gain Crossover. The phase margin at 
a frequency of gain crossover. 

Phase Margin, Minimum. In an absolutely stable feed- 
back loop, the least value of phase margin for those fre- 
quencies at which loop gain is between prescribed 
values. 

Position Constant. The reciprocal of the position error 
coefficient. 

Position Error Coefficient. The constant term of the 
power series expansion of the difference transfer ratio. 
Time Response. The output of an element or system as 
a function of time in response to a specified input. The 
state of the system at the time of application of the in- 
put must be specified. 

Velocity Constant. The reciprocal of the velocity error 
coefficient. 


Velocity Error Coefficient. The coefficient of the linear 
term in the power series expansion of the difference 
transfer ratio. 
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PGAC News 


ADMINISTRATIVE COMMITTEE 
MEETING 


The pertinent sections of the Minutes of 
the Administrative Committee Meeting, 
held on August 20, 1957 in the Fairmont 
Hotel, San Francisco, Calif., are given here 
in the belief that PGAC members should 
know what the Group policies and activities 
are and how they are being formed. Ques- 
tions and suggestions for future Group ac- 
tivities are always welcome. (The Adminis- 
trative Committee members are listed on the 
inside of the front cover.) 


1) Gene Grabbe reported on the organiza- 
tion of the PGAC. 


The program for the coming IRE Na- 
tional Convention in New York was dis- 
cussed. Gene Grabbe indicated that he 
would like Dave Lindorff of the University 
of Connecticut to serve again as program 
chairman for the Convention since he did an 
excellent job last year. It was suggested that 
he be aided in this work by a representative 
in the West, and Tom Stout accepted this 
assignment. Larry Cumming suggested that 
Dave Lindorff get further aid and Tom 
Mahoney accepted the assignment in the 
Boston area. It was agreed that Vic Azga- 
petian be asked to do the same in the New 
York area. It was felt that these four would 
compose an excellent program for this very 
important convention. The mailing list that 
has accumulated from previous arrange- 
ments of PGAC programs, and was used by 
Art Hopkin in arranging the WESCON 
sessions, will be sent to Dave Lindorff for his 
use. 

The committee assignments were then 
enumerated by Gene Grabbe: Awards— 
George Biernson; Meetings—Dave Lindorff ; 
1957 WESCON—Art Hopkin; 1958 Na- 
tional Convention—Dave Lindorff; Chap- 
ters—Bob Wilcox; Papers—George Axelby; 
1958 WESCON—to be assigned; Newsletter 
—John Ward; Publicity—to be assigned; 
1959 PGAC National Conyvention—Louis 
Wadel. 

The National Convention planned for 
the fall of 1959 was discussed to some extent. 
The desirability of holding it in the Dallas 
area and joining it to the Simulation Con- 
ference or the Southwestern IRE Conference 
for Region 7 was considered, but no definite 
decision was reached. Louis Wadel accepted 
the responsibility of surveying the situation, 
and submitting a proposal or alternate pro- 
posals. These proposals should be made 
available to members of the PGAC Adminis- 
trative Committee in time to be reviewed 
before the next meeting of the PGAC Ad- 
ministrative Committee during the New 
York Convention. 


2) Reports by Committees and Officers. 


Papers Committee—George Axelby de- 
scribed the mechanism of reviewing papers. 
He acts as a coarse filter and rejects obvi- 
ously improper material. The rest he sub- 
mits to three reviewers for comments. In 
most cases papers are returned to the author 
for modification. All this illustrates the time 
lag between receipt of manuscript and pub- 


lication. Louis Wadel suggested D. J. Sim- 
mons of Convair, Fort Worth, Tex., as re- 
viewer. 

Although advertising in the TRANSAC- 
TIONS is an approved IRE policy (at double 
the cost of PROCEEDINGS advertising), the 
committee felt that PGAC should not solicit 
for advertising at this time. However, it was 
agreed that Institute listings should be pro- 
moted. A list of organizations will be pre- 
pared by George Axelby. John Salzer sug- 
gested that a copy of the TRANSACTIONS be 
sent along with the requests for listing. A 
letter will be prepared by George Axelby 
and sent to Gene Grabbe and Larry Cum- 
ming for approval before being mailed out 
to the various organizations. 

Newsletter Committee—George Axelby 
was asked to transmit the information to 
John Ward concerning this topic. It was 
agreed that a Newsletter is-a very worth- 
while project for promoting the spirit and 
purposes of PGAC. Numerous suggestions 
were made as to what the Newsletter should 
include. It should include actions of the ad- 
ministrative committee, previews of the 
coming TRANSACTIONS, chapter news, per- 
sonal news items, etc. The Newsletter should 
be brief and well-written. 

Secretary-Treasurer Report—John Salzer 
reviewed the financial and membership 
status of PGAC. The organization has 
$9,319.91 which seems an excellent state of 
affairs considering the brief history of 
PGAC. Membership now stands at a total 
of 2568 of which only 75 are unpaid. Gene 
Grabbe made a correction of this figure to 
2570 by transferring to John Salzer $4 he 
had collected from newly recruited PGAC 
members. 

Chapters Committee—With Bob Wilcox 
absent, committee members reported on the 
various chapters. 

The formation of the San Francisco or 
Bay Area Chapters were affected by L. K. 
Lee leaving this area. Fortunately, Art 
Hopkin of the University of California 
agreed to initiate work on the San Francisco 
Chapter. George Axelby reported on the 
Baltimore Chapter which had two combina- 
tion dinner and technical meetings last year 
attended by about 65 and 45 persons, re- 
spectively. This chapter plans five meetings 
in the next season, some of which will be co- 
ordinated with the local AIEE Technical 
Group on Feedback Control. Tom Mahoney 
reported on the Boston Chapter which had 
five meetings last year with the attendance 
ranging from 35 to 125. He thought the suc- 
cess of the meetings was closely related to 
the prominence of the speaker. 

Louis Wadel reported on the Dallas-Fort 
Worth Chapter. He thought the chapter’s 
attendance had been somewhat low, but had 
great hopes for future interest and growth. 
The present chapter chairman is D. J. 
Simmons. 

John Salzer reported on the Los Angeles 
Chapter of which he had been chairman to 
this time. Meeting attendance ranged froma 
dozen to over 300. Several of the meetings 
had been with other groups and were pur- 
posefully held at various areas of the Los 
Angeles Section. The chapter has acquired a 
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real momentum and the coming year is ex- 
pected to be even more successful. 

It was noted that a new chapter was 
formed in Philadelphia. 

Consideration was given to the lack of a 
PGAC chapter in the Detroit area. John 
Salzer indicated that he might write a letter 
to A. C. Hall, Director of Bendix Research 
in Detroit, concerning chapter formation. 

The Affiliate Plan—The IRE Affiliate 
plan was discussed. Larry Cumming noted 
that the mechanism for affiliate membership 
in the IRE is completely set up. The plan 
has not gained momentum yet and there are 
actually very few affiliate members so far. 


3) Meetings. 


The Nonlinear Symposium—This sym- 
posium, held in San Francisco prior to the 
committee meeting, attracted a registration 
of 114 persons plus speakers and panel mem- 
bers. The luncheon was attended by more 
than 100 persons. Cost of the convention 
was $75 for the meeting room and $573.75 
for program printing. 

AIEE October, 1957, Conference on Com- 
puters in Control, Atlantic City—Gene Grabbe 
reported that PGAC was a participating 
organization and contributed about a half 
dozen papers. 

ASME, April, 1958, Conference on Opti- 
mizing Control, Delaware and AIEE, April, 
1958, Conference on Automatic Techniques, 
Detroit—These are two meetings which 
PGAC is interested in supporting. 

The National Control Conference, Oc- 
tober, 1959—This conference has been dis- 
cussed above and will be a PGAC-sponsored 
affair. Louis Wadel is in charge of plans. 

IRE National Convention, March, 1958— 
This convention has also been discussed 
above. Dave Lindorff is in charge of the pro- 
gram. 


4) Other Business. 


Gene Grabbe announced the formation of 
a National Automatic Controls Committee 
to represent the U. S. in the International 
Federation of Automatic Controls. NACC 
is a vehicle to achieve cooperation among 
five American groups in the automatic con- 
trol field: IRE, AIEE, ASME, AICHE, and 
ISA. The next meeting of the IFAC was 
held in Paris on September 10, 1957 with the 
following representation for NACC: J. 
Lozier for IRE, Harold Chestnut for AIEE, 
and Rufus Oldenburger for ASME. 

To provide NACC with general working 
funds, the PGAC administrative committee 
voted a contribution of $300 which was 
matched by other member societies of 
INNCG: 

Mr. Petrie, of the IRE Student Relations 
Committee, brought up the question of get- 
ting students interested in automatic con- 
trols at the high-school level. He suggested 
that at a convenient time PGAC might con- 
sider appointing a committee to work out 
plans for achieving this purpose. One such 
plan might be to set up instructions im- 
provising some simple controls that a stu- 
dent could build at school or in his home. 
Another approach might be the provision for 
certain competitive tasks with prizes. 
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Region 1 


Binghamton 


Bernstein, R. 
Bosman, E, H. 
Dibb, G. 

Evans, B. O. 
Gibson, R. G. 
Hamburgen, A. 
Hemstreet, H. S. 
Ivy, R. C. 
Johnson, G. W. 
Kercher, R. B. 
Kilmer, F. G. 
Kovalchick, N. 
Lavender, R. W. 
Lohman, I. H., Jr. 
Shatz, J. R. 
Sitterlee, L. J. 
Van Horn, J. F. 


Boston 


Ackerman, S. 
Aiken, R. T. 
Alden, J. M. 
Alexander, C. S., Jr. 
Allen, J. 
Anderegg, J. S. 
Annis, E. W., Jr. 
Antul, J. J. 
Applegate, C. E. 
Apsokardos, E 
Barry, J. G. 
Batchelder, L 
Beaudette, C. G 
Beaumariage, D. C. 
Beecher, A. E. 
Bell, C. G. 
Benkley, F. G. 
Bennett, H. W. 
Bennett, R. K. 
Biernson, G. A. 
Blanchard, R. L. 
Bliss, Z. R. 
Bloomfield, B. 
Bosch, F. M. 
Bradley, C. L. 
Brew, R. A. 
Brindley, W. E. 
Brooks, P. R., Jr. 
Brown, B. D. 
Brown, G. S. 
Brown, J. B. 
Bullard, A. H., Jr. 
Burwen, R. S. 
Buxton, E. T., Jr. 
Capen, E. B. 
Carrier, G. T. 
Chin, B. M. 
Chul. He 
Clapp, C. W. 
Clements, D. F. 
Connelly, M. E. 
Cook, C. G. 

Cox; G. C. 
Dandreta, W. 
Davidow, W. 
Demrow, R. I. 
De Russo, P. M. 
Dickson, A. W. 


Dratch, J. E. 
Dworshak, F 
Edie, J. L. 


_ Engel, J. S. 
Bie J. K.- 


I 


Estes, S. E. 

Evans, R. R. 

Fallows, E. M. 

Farmer, J. W. 
TEES Eley yam Von be 
Farrah, H. R. 

Fegley, K ie 
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Roster of PGAC Members 


Listed by IRE Regions and Sections as of January 15, 1958 


A year ago, the membership of the PGAC was published in the February, 1957 issue 
of these TRANSACTIONS. Since that time the membership has more than doubled. To 
indicate the growth of the PGAC and the engineers who identify themselves profession- 
ally with automatic control, the membership directory is again published by regions and 


section including members of foreign countries. 
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Fertig, K. 


Fitzmorris, M. J., Jr. 


Frank, W. I. 
Freeman, L. D. 
Freudberg, R. L. 
Fricks, R. E. 
Galagan, S. 
Gallagher, E, F. 
Garber, T. A. 
Gitelman, E. 
Goldberg, D. 
Gordon, B. M. 
Goulder, M. E. 
Grossimon, H. P. 
Guethlen, V. J. 
Gutwill, L. A. 
Hamza, M. H. 
Hansen, J. T. 
Heaviside, M. G. 
Hersh, A, I. 
Heuchling, T. P. 
Hillman, H. D. 
Hills, F. B. 

Hills, W. L. 
Holubow, F. E. 
Hoogasian, L. 
Hopkins, A. L., Jr. 
Howard, R. A. 
Huffman, D. A. 
Hulburt, H. J., Jr. 
Iffand, J. J. 
Jamgochian, E. 
Johnson, E. P., Jr. 
Kaiser, J. F. 
Kaplan, M. L. 
Karmel, P. R. 
Kasvand, T. 
Kaye, M. G. 
Kellner, R. 
Ket CL: 
Khambaty, M. B. 
Kilroy, H. P. 
Kleinrock, L. 
Kodis, R. D. 
Korn, S. J. 
Kotliar, A. 
Kremheller, W. G. 
Krulee, R. L. 
Lechner, R. J. 
Leonard, C. E. 
Leonard, R. R. 
Levy, D. M. 
Lincoln, A. J. 
Lindsay, G. A. 
Lovell, R. W. 


Mahoney, T. F. 
Marino, A. S. 
Mark, R. B. 


McMuttrie, Dies 
McNamara, J. N. 
Meditch, J. S 
Melanson, F. J. 
Mercer, W. R 
Merrill, B. 
Miles, R. J. 
Milias, E. 

Miller, H. A. 
Minsky, M. L. 
Misek, V. A. 
Morgenstern, sae 
Morris, R. V. 
Mortensen, R. E. 
Murano, L. 


Naylor, 4 Ks. 
~ NY ae 


vy. EAA 


Newton, G. C., Jr. 
Nielsen, C. E., Jr. 
O'Brien, D. G. 
Oettinger, A. G. 
Olsson, E. A. 
Orenberg, A. 
Osman, M.S. 
Packard, R. H. 
Palmer, P. J. 
Pantazelos, P. G. 
Parker, D. A. 
Pastan, H. L. 
Pease, W. M. 
Perron, R. R. 
Perry, B. D. 
Perry, K. E. 
Petteruti, A. J. 
Phillips, D. R. 
Platt tig. 
Poor, V. D. 
Prival, H. G. 
Pugh, Ae hg Ill 
Pughe, E. W., Jr. 
Reichard, R Ww. 
Repella, N. 
Rittenburg, S. E. 
Roch, M. E. 
Roediger, F. E. 
Rojak, F. A. 
Rona, T. P. 
Rowell, W. G. 


Scheidenhelm, R. 
Schochet, J. R. 
Schorr-Kon, J. J 


Schramm, M. W., Jt 


Seifert, W. W. 
Shansky, eS 
Sheingold, D 

Shortell, A. V. rae 
Shulman, L. 

Sims, B. S. 
Sinclair, D. B. 
Small, R. M. 
Smith, B. K. 
Smith, R. L. 
Smith, T. B. 
Snyder, D. C. 
Soderstrom, R. E. 
Solomonoff, R, J. 
Spink, P. G. 
Steinberg, C. A. 
Strassberg, D. D. 
Susskind, A. K. 
Swonger, C. W. 
Taylor, H. P. 
Teixeira, N. A. 
Thayer, W. D. 
Lsa0, C.K. . 
Tunnicliffe, W. W. 
Vacca, R. H. 
Valpey, R. S., Jr. 
Vander Velde, W. TE. 
Ventura, J. J. 
Viggelaar, H. V. D. 
Vincent, G. D. 
Walsh, es Pr 
Ward, J. E. 

Wexler, H. 4k, 

el ae he WA. Jr. 
White, D. 
Whittaker, H. ~ jr, 


“Wilcox, R. B. 


Wilkie, L. E. 
Williams, S. B. 
Wing, J 

Wiren, J. 
Woodman, R.¥F. 
Woodruff, T. E. 
Woodward, J. H 
Young, F. M. 
Zieman, H. E 


Buffalo-Niagara 


Aines, F. G. 
Archibald, W. R. 
Beneke, J. 
Hrank i. B. 
Gray, R. B. 
Green, V. H. 
Grose, C. W. 
Halsted, G. P. 
Hawes, EaG: 
Hayman, R. A. 
Johnson, A. B. 
Meister, L. E. 
Michaelis, T. D. 
Newton, D. J., Jr: 
Novo, R. B. 
Powell, F. D. 
Rochin, D. W. 
Ruda, E. V. 
Savi, L. 
Walbesser, W. J. 


Connecticut 


Bourret, C. J. 
Brockett, R. I. 
Brooks, F. A., Jr. 
Brownell, R. M. 
Curtiss, R: H. 
Delany, E. B. 
Dellenbaugh, F. S. 
Dolan, O. F. 
Farwell, R. P. 
Flynn, ©. F. 
Georgi, E. A., Jr. 
Gilchrist, E. S. 
Gouyet, J. P. J. 


Gray, N. L. 
Gutterman, A. Z. 
Haas, V. B., Jr. 

Hall, B. A. 
Harrington, C. F., Jr. 
Henry, H. E. 
Hubbard, R. G. 
Jensen, H. T. 
Kidwell, N. E. 


Kosowsky, L. H. 
Krauss, H. G. 
Lamb, J. J. 
Leuze, W.F. 
Lindorff, D. P. 
Lutz, C. H. 
Maltais, W. E. 
Mapes, T. J. 
Miller, E. F. 


Montgomery, R.YH., Jr. 


Murphy, R. E., Jr. 
Murray, G. J., Jr. 
Narendra, K. S. 
Nelson, G. E. 
Northrop, R. B. 
Norton, C. A. 
Orell, F. M. 
Palmero, A. 
Passera, A. L. 
Penney, A. W. 
Petruchelly, V. J. 
Plehaty, S. L. 
Poland, W. L. 
Sherman, P.M. 
Sontheimer, C. G. 
Stephanz, G. H. 
Stockman, W. E. 
Stoker, W. C. 
Sunega, E. M. 
Sweet, R. D. 
Swift, W. C. G. 
Torrance, Ap dal 
Townsend, F. H. 
Williams, i But. 
Zweig, F. 


Erie 
Lempert, J. 


Ithaca 


Alcaide, H. D. 
Criswell, T. 
Dockwiller, W. R. 
Hufnagel, R. E. 
Jackson, A. S. 
Malone, D. M. 
Martin, A. R. 
Mayer, H. F. 
McLange, T. 
Meserve, W. E. 
Phillips, TD. Az 
Scudder, H. J., III 
Modde) nts 


Rochester 


Berch, W. H. 
Brown, G. A. 
Chesna, J. 
Dutcher, B. C. 
Bilis, LT. E., Jr. 
Enslein, K. 
Faust, A. C. 
Federici, J. T. 
Formicola, A. F. 
Haskell, J. W. 
Heit, J. C 
Johnson, R. M. 
Merle, C. W. 
Morse, J. E. 
Rodatus, E. J. 
Ruppert, J. R. 
Shalloway, A. M. 


Sheehan, J. F. 
Shepard, W. 4H. 
Stone, D. J., Jr. 
Thalner, R. R. 
Traub, R. A. 
Trott, M. 


Rome-Utica 


Baldridge, B. H. 
Glaser, G. J. 
Hatfield, J. P. 
Kabrisky, M. J. 
Leney, H. G. 
Rossy Pai: 


Schenectady 


Acton, E. S. 
Borner, E. F. 
Buchhold, T. A. 
Buscher, R. G. 
Chestnut, H. 


Fanuele, F. J. 
Kirchmayer, L. K. 
Liebowitz, B. H. 
Lippitt, D. L. 
Pester, R. F. 
Rothe, F. S. 
Shuey, R. L. 
Stromer, P. R. 
Stutt, C. A. 
Wright, W. H., Jr. 


Syracuse 


Bessette, D. U. 
Brady, D. J. 
Buchta, J. C. 
Cottle, D. W. 
Cowan, F. C. 
Edwards, K. A. 
Eibeck, A. C 
Fehlau, C. E 
Grisetti, R. S. 
Heartz, R. oh 


AS ureller, J i 
McCarthy, J. J. 


Neelands, L. J. 
Russell, J. B., Jr. 


: 
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Shuart, O. H. 
Smith, G. H. 
Stabler, E. P. 
Vaughan, J. A. 
Zuroff, C. A. 


Western Massachusetts 


Anderson, J. M. 
Lovell, B. W. 
Luoma, R. A. 
Mahar, T. F. 
Patch, R: J. 


Region 2 


Long Island 


Adise, H. H. 
Agree, I. 

Amato, J. A. 
Austin, G, A. 
Barber, E 
Bargeski, A. J. 
Beck, D. H. 
Behn, E. R. 
Bennett, H. A. 
Boxer, R. 
Buehrle, W. E., Jr. 
Burgess, E. G., Jr. 
Burr, R. P. 


Cardwell, ROA 

Caruthers, F. P. 
Chapman, P. W. 
Chartoff, P. 


Comenzo, L. F. 
Corrado, V. M. 
Costa, T. A. 
Crosby, M. G. 
Dettinger, D. 
Detwiler, S. P. 
Di Toro, M. J. 
Dmytrasz, J. N. 
Doersam, C. H., Jr. 
Dyer, J. N. 
Egnuss, E. M. 
Ellis, P. H. 
Engelson, H. R. 
Feay, D. I. 
Firth, L. G., Jr. 
Fishbein, M. 
Fonseca, A. P. 
Forman, L. E. 
Francheschini, J. B. 
Frank, F. P, E. 
Freed, A. 
Freeman, H. 
Fried, G. 
Friedman, E. D. 
Fromm, W. E. 
Fullerton, W. T. 
Gennard, R. D. 
Gilbert, R. C. 
Giordano, S. 
Glixon, H. R. 
Glowacki, J. 
Gneses, M. I. 
Gold, M. R. 
Goodstein, J. 
Gordon, R. L, 
Gorelick, G. 
Grabbe, D. 


Klein, R. 
Knocklein, H. P. 
Knox, R. W. 
Kokkinos, C. 

Kaird, J. A., III 
Lampathakis, K. E. 
Laspina, C. A. 
Laudon, H. C. 
Lawrence, C. W. 

Lee, C. H. 
Lemanczyk, J. C. 
Lenefsky, S. 
Levenstein, H. 
Levinson, E, 
Lieberman, A. G. 
Marcinkowski, E. R. 
Marston, R. S. 
Match, M. J. 
McPherson, D. L., III 
Meirowitz, R. L. 
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Menes, J. 
Mesecher, M. G. 
Moritz, F. G. 
Murphy, G. N. 
Murphy, R. B. 
Murtagh, J. B. 
Odessey, P. H. 
Osder, S. S. 
Osterndorf, J. F. 
Otten, H. C 
Packer, L. 
Padovano, R. J. 
Papke, F. E. 
Pearsall, C. H., Jr. 
Perliss, R. E. 
Peterson, iBE, Gy 
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Whaley, W. W. 
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Silvera, R. V. 
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Aaron, M. R. 
Abraham, R. P. 
Acker, J. L. 
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Anderson, N. E. 
Antonazzi, F. J. 
Archidiancono, J. 
Bahls, W. E. 
Bailey, E. M. 
Baumann, A. H. 
Bearman, A. L. 
Benz, B. D. 
Bogner, I. 
Boubli, E. J. 
Bowker, M. W. 
Brendle, T. A. 
Brown, A. T., III 
Brown, C. S. 
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Mathews, M. V. 
Mayo, J. S. 
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McCrory, J. R 
McIntyre, J. W. 
Meinholtz, H. J. 
Metzger, J. E. 
Morrow, S. R. 
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Mueller, P. L. 
Mulligan, J. H., Jr. 
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Ossanna, J. F., Jr: 
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Myers, F. G. 
Obermayer, R. W. 
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Parker, R. L.,. Jr. 
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Cherry, W. H., Jr. 
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Littleton, W. W. 
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Stephens, T. L. 
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Riggs, D. L. 
Rogers, R. F. 
Rogers, R. S. 
Roop, R. W. 
Rudnick, J. J. 
Rudofsky, S. 
Rydz, J. S. 
Seawell, W. N. 
Selinsky, J. J. 
Sevian, E. A. 
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Solomon, S. 
Sommer, W. G. 
Sorkin, C. S. 
Steudel, W. R. 


Sunstein, D. E. 
Sweeney, C. W. 
Tonooka, S. 
Tweet, B. O., Jr. 
Walker, H.R. 
Weiner, J. R. 
Weisenberger, A. J. 
Weiss, E. 

Wentz, J. L. 
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Willis, J. T. 

Wills, W. P. 
Wirtz, E. L. 
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Wolfson, H. S. 
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Wood, R. F., Jr. 
Wooten, L. B. 
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Ziegler, I. A. 
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Andrews, R. E. 
Black, J. H. 
Branscom, G. A. 
Burlingame, C. W. 
Cockrell, W. D. 
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Dial, E. W. 
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Polak, H. 
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Flowers, H. L. 
Haas, D. L. 
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Hurley, W. A. 
Lambert, C. O. 
Meilander, W. C. 
Miller, J. H. 
Nelson, R. L. 
Penniman, I. B. 
Ryburn, P. W. 
Stahl, M. D. 
Toman, W. J. V. 
Yarosh, N. P. 
Yochelson, S. 
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Bennett, P. E. 
Harvey, H. B. 
Hoechner, I. L. 


Knausenberger, G. E. 


Lawther, J. M. 
Oblinger, J. T. 
Seeley, R. M., Jr. 
Wolfe, R. E. 
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Berg, D. F. 
Blasberg, L. A, 
Brown, D. L. 
Butterworth, G. S. 
Colclaser, R. A. 
Dale, W. L 
Doerr, W. H. 
Engelmann, R. H. 
Forbes, H. R. 
Fortier, R. E. 
Georger, L. J. 
Herrin, C, B. 
Hocking, L. J. 
Keene, L. C. 
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Kinkaid, J. C. 
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Miller, H. J. 
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Phillips, W. E., Jr. 
Saltzer, C 
Shepherd, B. R. 
Tame, J.S 
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Bain, Tf. D:, Jr 
Burgener, R. C. 
Chope, H. R. 
Cohen, D. 
Conlon, R. J. 
Cosgriff, R. L. 
Cummins, M. M. 
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Peterson, L. S. 
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Thompson, J. P. 
Thompson, N. P. 
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Wichmann, Tae 
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Detroit 


Amber, P. S. 
Barcus, R. 
Blackwell, W. A. 
Bozoian, M. 
Brown, L. R 
Brunais, E. G. 
Bublitz, A. T. 
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Johnson, R. R. 
Kazda, L. F. 
Klem, R. F. 
Koziara, E. H. 
Lazarus, D. H. 
Lindahl, C. E. 
McGlinn, E. J. 
McGregor, A. D. 
Morgan, B. S., Jr. 
Mort, V. A. 
Nakagawa, N. 
Nixon, J. D. 
Olson, R. G. 
O’Neal, R. D. 
Patton, H. W. 
Raveh, Db. E. 
Rehwoldt, T. V. 
Retko, E. 
Robinson, D. F. 
Roth, P. F. 
Scott, D. L. E. 
Seleno, A. A. 
Sims, R. C, 
Smith, B. A. 
Smith, D. F. 
Smith, W. 
Strand, J. 
Sutton, W. A. 
Taplin, L. B. 
Thompson, R. H. 
Tokad, Y. 
Tubbs, R. J. 
Turkish, M. C. 
Webber, R. C. 
White, O. 
Willett, R. M. 


Pittsburgh 


Aronson, M. H. 
Baer, W. E. 
Bhavnani, K. H. 
Blewitt, D. D. 
Boyd, J pel 
Caywood, W. P., Jr. 


Mathias, R. 
Moss, A 
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Murray, D. W. 
O'Donnell, J. J 
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Schwindt, A. J. 
Spriggs, L. A., Us 
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Strull, G. 
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Butler, J. A. 
Ewing, D. J., Jr. 
Fuller, L. E. 
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Murley, E. M., Jr. 
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Alessio, S. A. 
Arndt, L. K. 
Auth, L. V., Jr. 
Axelrod, L. R. 
Barnett, W. T. 
Bielenin, A. 
Bold, N. a 
Boyd, D. M., Jr. 
Brietze, F. 


Bullen, C. V. 
Burtness, R. W. 
Carter, R. 
Cermak, C. W. 
Chorney, P. L. 
Chulsky, I. 
Condron, W. Be 
Cruze, Be 
Davoust, D. i 
Dennis, R. J. 
Deterville, R. J. 
Dikinis, D. V. 
Druz, W. S. 
Ebstein, B. 
Eilers, C. G. 
Fenwick, W. D. 
Ferre, G. E. 
Foster, G. E. 


Gehrke, W. C. 
Gerlach, A. B. 
Glyptis, N. 
Gozner, S. M., Jr. 
Greenberg, C. J. 
Gregory, E. C. 
Heid, K. K. W. 
Isolampi, G. E. 
Jenness, R. R. 
Kortright, F. U. 
Kott, W. O 
Kreer, J. B. 
Kuchenbecker, R. A. 
Lafferty, V. C. 
Lesnick, R. 


Martin, is W., Jr. 
Masteller, Bela 
Mercure, R. 
Merrifield, L. A. 
Meyer, A 
Mittelmann, E. 
Motch, J. F. 
Newman, J. B., Jr. 
Noges, E. 

Olson, J. C. 
Paterson, W. L. 
Rebman, W. H. 
Rockwood, C. C. 
Ruther, H. E. 
Sayles, H. L. 
Severns, R. A. 
Shewan, W. 
Sigborn, T. W. 
Slana, M. F. 
Sommeria, M. R. 
Thielen, L. R. 
Thomas, R. G. 
Van Bosse, J. G. 
Vandling, G. C. 
Van Ness, J. E. 


Van Valkenburg, M. E. 


Veazie, W 
Venne, R, E. 
Veras, J. A. 
Verbanec, W. R. 
Vitous, J. P. 
Warshawsky, J. 
Waters, J. I. 
Wavering, A. J. 
Weissman, R. M. 
White, E. S. 
Woolums, L. L. 
Ye, J. W. F. 


Fort Wayne 


Brady, F. H. 
Cronkwright, R. E. 
Emery, R. C 

Fife, D. W. 
Johnson, ae we 
Kaplan, R. 
Klingler, E. H. 
Mason, C. F. 
Nelson, G. L. 
Richeson, W. E. 
Solomon, R. M. 
Williams, W. J., Jr. 


Indianapolis 


Albon, R., Jr. 
Braunagel, M. V. 
Cannon, C. D. 
Carpenter, M. M., Jr. 
Cross, K. R. 
Ennis, F. L. 
Evans, R. A. 
Gibson, J. E. 
Groves, C. R. 
Halsey, M. R. 
Israel, J. D. 
Johnson, T. L. 
Jones, M. E. 
Luisi, J. A. 
Marshall, H. B. 
Ogborn, L. L. 
Skinner, G. M. 
Spencer, J. L. 
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Tous J. 
Whipple, R. L. 
White, S. A. 


Louisville 


Ebaugh, D. P. 
Kain, R. Y. 
Kwo, T. T. 


Milwaukee 


Arakelian, G. P. 
Bennett, T. H. 
Camps, T. F. 
Carlson, A. W. 
Diggelmann, H. 
Frankos, D. T 
Gessner, U. 
Graham, J. D. 
Jensen, K. S. 
Kerske, J. F. 
Laurents, V. T. 
Leary, R. W. 
Limpel, E. J. 
ind? sR. 
Lindemann, A. W. 
Lofy, F. J. 
Louzader, J. C. 
Mayer, H. J. 
Mezger, J. P. 
Morin, D. C.., Jr. 
Neider, F. V. 
Norum, V. D. 
Pierce, R. L. 
Pierick, K. R. 
Rekoff, M. G., Jr. 
Smithy CoC. 
Stephan, R. R. 
Weissenborn, H. E. 
West, B. M. 


Omaha-Lincoln 


Bashara, N. M. 
Terry, R. W. 


South Bend-Mishawaka 


Burke, @. A. 
Bymbaent RS ee Bh 
Chen, G. K 

Cooney, J. i ' 
Hansen, A. G., Jr. 
Hoffman, C. H. 
Minnick, V. P. 


Twin Cities 
Adams, G. E. 
Alderson, R. C. 
Alfsen, G. F. 
Alterman, F. J. 
Balzart, EB: Js Jr: 
Bartlett, V. W. 
Benassi, D. A. 


Carlson, R. A. 
Evans, J. M. 

Gise, F. G., Jr. 
Gustafson, insta WN 
Hardenbergh, G. A. 
Harvey, P. O. 
Hird, F. S. 
Hulstrand, B. E. 
Inman, T. F. 
Kershaw, J. A. 
Kiene, R. C. 
Knoblauch, A., Jr. 
Kukuk, H. S. 
Lahue, P. M. 
Lanzkron, R. W. 
Lease, L. R. 

Lode, T. 

Ludwig, J. T. 
Luik, I. 
Markusen, D. L. 
Maze, R. O. 
McLane, R. C. 
Moe, W. J. 
Murray, R. L. 
Nellis, W. M. 
Noble, D. S. 
Nordstrom, J. E. 
Ormsby, R. D. 
Pierce, A. J. 
Ramamoorthy, C. V. 
Reed, M. W. 
Rowland, C. A., Jr. 
Schellenberg, A. G. 
Schiebe, F. R, 
Schuck, O. H. 
Senstad, P. D. 
Smith, D. A. 
Stone, N. T. 
Storm, J. F. 
Swanholm, W. J. 
Swanlund, G, D. 
Vogel, J. P. 
Windsor, R. N. 
Woll, L. J. 


Region 6 


Alamogordo-Holloman 
Carlson, H. A. 
Chandler, D. P. 
Koellner, K. H. 


Liston, D. H. 
Pykosz, T. F. 


Dallas 


Anderson, R. P. 
Askew, W. J., Jr. 
Bailey, R. W. 
Barnett, M. L. 
Braun, C. 

Buehrle, C. D. 
Dubose, G. P. 
Gentry, J. A. 
Humke, F. O., Jr. 
McDonald, M. 
Miller, N. D. 
Osnaschek, F. J. 
Pittman, P. D., Jr. 
Scammel, B. C. 
Sparkman, W. A., Jr. 
Stanton, A. N 
Tatum, F. W. 
Teasdale, A. R., Jr. 
Wadel, L. B. 
Wilhelm, E. S. 
Ziemer, D. R. 


Denver 


Bergman, R. R. 
Capehart, M. E. 
Cook, E. E. 
Daniels, W. H. 
Hart, W. G., Jr. 
Howell, R. G. 
Messler, F. J. 
Mielziner, W. 
Morroni, D. J. 
Ostwald, L. T. 
Spafford, L. J. 
Stacey, D. S. 


El Paso 


Bullock, R. H. 
Martin, T. A. 
Moore, T. E. 
Moss, E. C. 
Rojas, A. M. H. 


Fort Worth 


Delaney, J. S. 
Evans, W. L. 
Forester, J. R. 
Hines, R. L. 
Jiles, C. W. 
Jones, L. R. 
Perry, E. R. 
Seale, I. A. 
Simmons, D. J. 
Slagle, G. M., Jr. 
Watkins, O. E. 


Houston 
Bueys Jake, ts 
ads, iR&) ire 
Easterling, M. F. 
Frobese, C. W. 
Groendycke, A. R. 
Hargett, J. T. 
Keating, L. M. 
Keeling, W. E. 
Kuykendall, W. E. 
Lockwood, R. M. 
Navarro, S. O. 
Sherrill, B. 
Tasini, B. 
Terao, M. 


Kansas City 


Annis, J. W. 
Brown, G. M. 
@lum, Ll. E: 
Halijak, C. A. 
Lobb, C. W. 
Rauch, E. L. 
Reiser, R. R. 
Roehm, R. A. 
Simonds, R. L. 
Stout, H. L. 
Wilcox, J. V. 


Little Rock 
Anthony, C. A. 
Bodie, M. W. 
Butler, B. L. 
Cannon, W. W. 


Lubbock 
Perkins, C. S; 


New Orleans 


Cronvich, J, A. 
Dietz, W. R. 
Drake, R. L. 
Hodge, B 

Maus, L. C. 
Smith, L. G. 
Weiser, C. C., Jr. 
Welsh, M. J., Jr. 


Oklahoma City 


Glennon, R. E. 
Grubbs, C. E. 
Ledbetter, R. P. 
Vlay, G. J. 


St. Louis 
Abkemeier, L. R. 


Barr, R. K. 
Cassidy, R. E. 
Dawson, R. T. 
Ebert, W. J., Jr. 
Fiedler, G. J. 
Hamilton, S. 
Herchenreder, H. 
Hudgens, L. L. 
Lago, 

Malsbary, ie S. 
Mayer, M. F., Jr. 
McAninch, CG. igh 
Min, H. S. 
Mohrman, R. F. 
Mutchek, J. H. 
Norman, C. F. 
Palmer, J. A. 
Reed, D. L. 
Sayer, J. D. 
Sceherz,|Conic 
Sheehan, J. S. 
Sudfeld, C. C. 
Twombly, J. W., Jr. 
Weaver, M. D. 
Zaborszky, J. 


San Antonio 


Bostick, F. X. 
Fisher, H. V., Jr. 
Mirsch Gs 0% 
Hoffman, A. A. J. 
Phillips, J. P. 
Reinhard, E, A. 
Saadeh, D. 


Shreveport 


Averre, C. W. 
Gordon, E. 
Long, F. V. 


Tulsa 


Brashear, R. T. 
Day, C. E. 
Freeman, L. R. 
Fox, D.N. 
Greening, J.P. 
Labarthe, L. C. 
Peterson, G. 
Piety, R. G. 
Rowley, R. G., Jr. 
Silverman, D. 
Sykora, G. E. 
Wolkov, D. 


Region 7 


Albuquerque-Los Alamos 


Baker, D. A. 
Carlson, W. F. 
Foster, F. N. 
Helgeson, B. P. 
Buy. 0: 
Koschmann, A. H. 
Koskela, A. C. 
Pace, T. L. 
Powell, R. L. 
Ray, H. K. 
ostin Dai. 


China Lake 
Creusere, M. C. 


Poulson, W. A. 
Starr, E. A. 


Los Angeles 


Adomian, G. 
Akers, R. C. 
Akin, P. A. 
Albrecht, A. 
Alfrey, P. R. 
Allen, D. H. 
Alper, S. M. 
Ambrose, J. R. 
Amstutz, M. F. 
Anderson, F. 
Anderson, M. J. 
Anderson, W. E. 
Andrews, L. V. A. 
Anzel, B. M. 
Arnold, he 
Aroyan, G. F. 
Arsenault, W. R. 
Aseltine, J. A. 
Auletti, R. L. 
Avrech, N. 
Baker, D. L. 
Baker, J. N. 
Barlett, F. R. 
Barnett, L. 
Barston, E. M, 
Bartholomew, H. R. 
Beard, G. W 
Beckwith, se W. 
Bekey, G. A. 
Belding, R. A. 
Bement, W. A. 
Bemis, PLO 1h, 
Berg, R.L., Jr: 
Berman, M. R, 
Bertorello, R. J. 


Bible, R. E. 
Bills, G. W. 
Bonney, R. B. 
Borgeson, P. W. 
Borsch, K. S. 
Bower, J. L. 
Brandt, R. 
Braun) BL. 
Braverman, D. J. 
Brewer, H. E. 
Broadwell, W. B. 
Bronstein, L. M. 
Brown, D. E. 
Bucher, F. X. 
Buland, R. N. 
Bunkowski, N. H. 
Burk, W. A. 
Burnsweig, J., Jr. 
Byck, D. M. 
Cahn, J. M. 
Callot, S. 
Campbell, G., Jr. 
Capps, J : 
Carlson, C. O. 
Chandaket, P. 
Chang, B. 
Chapman, C. W. 
Chernoff, D. P. 
Christensen, A. V. 
Cochran, E. D. 
Concus, P. 
Conway, J. P. 
Cooke, C. W., Jr. 
Corley, C. F. 
Cornell, J. R. 
Corvi, J. A. 
Cosgrave, S. J. 
Craven, W. A., Jr. 
Crowe, J. W. 
Curry, W.S.,- Jt. 
Cutting, E. 
Dains, H. O. 
Davis, H. G., Jr. 
Day, R. N. 
Deaux, F. J. 
Deen, J. R. 
Deming, A. F. 
Determan, J. D. 
Deuser, D. A. 

Di Bias, R. L. 
Dickinson, H. B. 
Diem, C. W. 
Diemer, F. P. 
Dietz, D. P. 

Di Giuseppe, A. 
Dinning, J. R. 
Doty, R. L. 
Drucker, A. N. 
Dzilvelis, A. A. 
Edelsohn, C. R. 
Eggeman, D. J. 
EHichwald, W. F. 


Eikelman, J. A., Jr. 


Ellis, D. O 
Engel, H. L. 
Eno, R. F. 
Epstein, S. 
Estrin, G. 
Evfimenko, A. 
Finley, W. A. 
Fischer, P. F. 


Forbath, F. Pe 
Fox, C. W. 
Foxman, E. 
Francis, T. F. 
Furumoto, N. 
Gabler, R. T. 
Garber, L. F. 
Garcia, E. M. 
Gardner, F. H. 
Gardner, F. M. 
Gauronskas, P. P. 
Gaylord, R. S. 
Gee, L. C. 


Glowalla, ie 
Goldstein, A. 
Gore, M. R. 
Grabbe, E. M. 
Graham, J. D. 
Green, R. D. 
Greensberg, E. L. 
Gronstrand, T, L. 
Gross, W. 
Gullatt, S..P., Jr. 
Gunning, W. F. 
Gustin, J. T. 
Hadden, F. A. 
Hailey, R. D. 
Hall Glas 
Halloran, R. G. 
Hanna, R. E. 
Harmon, W. G. 
Hassel, R. R. 
Hawkins, J 
Hayes, J. E. 
Helland, J. 
Heyliger, G. E. 
Hicks, A. R. 
Hinrichs, K. 
Hitchcock, R. W. 
Holmen, R. E. 
Hom, F, M. 


March 


1958 


Hoskinson, E. A. 


Hurley, R. F, 
Hutcheon, R. S. 
Izuel, A. (ee 
Jack, R. W. 
Jackson, C. F. 
Jackson, K. R. 
Jacobs, J. E. 
Jacobson, O. M. 
Janeway, R. K. 
Jenkins, L. E. 
Jeschke, A. W. 
Joerger, J. C. 
Johnson, M. G. 
Johnson, W. A. 
Jones, J. M. 
Kahn, B. S. 
Kaufman, F. H. 
Kaufman, S. 
Kawahata, B. I. 
Kazarian, D. G. 
Kennedy, C. J. 
Kennedy, F. D. 
Kennel, J. M. 
Keppel, R. A. 
Kern, W. W. 
Kerster, G. 


Kirsch, H. A. 
Kishi, F. H. 
Kitabayashi, T. 
Klimowski, F., Jr. 
Knox, R. V. 
Kobayashi, R. 
Krill, C. K. 

Kroy, W. H.., Jr. 
Kuerschner, H. 
Kukel, J. 

Kuntz, G. L. 
Lawrence, A. F., III 
Lebell, D. 


Lee, H. J. 
Leondes, C. T. 
Leone, W. C. 
Lessley, T. D. 
Levine, L. 
Levine, S. E. 
Levinson, R. M. 
Levy, E. C. 
Lewis, D. E. 


Lichterman, G. 
Lillibridge, E. H. 
Lindholm, C. R 
Lloyd, A. T. 
Lords, F. V. 
Louie, W. 
Lyon, D. B. 
Lyons, L. H. 
Makush, W. S. 
ager age Ao Re 
Manly, R. 
Manser, D. E., Jr. 
Margolis, M 
Maticich, J. R. 
Mayberry, L. A. 
Mayner, G. S. 
McCarthy, J. W. 
erence G. F. 
McCann, L. 
McGhee, R. B. 
McLarin, 
McLeod, M. Ge 
McRuer, D. T. 
McVey, I. M. 
Mehner, E. W. 
Metzner, H. E. 
Mileson, D. F. 
Miller, D. S. 
Mitsutmoi, T. 
Medlinski, W. M. 
Moise, N. L. 
Monroe, F. R. 
Moore, P. G. 
Morris, G. R. 
Morrison, A. I. 

__ Morton, W. B., Jr. 

- Mosher, W. W.., Jr. 

= Moss, J.F., Jr. 

_ Mundt, E. 

© Myers, W. 4r0 


Nelson, C. S., Jr. 
Nod 


Nuban, E. 
Nussbaum, O. N. 
Nuttall, H. V. 

4 O'Brien, | W.C. 

= Odlum, F. H., Jr. 
O’Hara, Cage 
Ohrenstein, S. B. 
- Olsen, J. D. 

- Olsen, L. ' 
Parker, A. 
= Parzl, R. C. 

; Patrusky, N. 
Perez, A. oO 
Perkins, L x M. 
Pernick, L. J. 


Peterson, E. J. 
Peterson, R. W. 
Phister, M., Jr. 
Poirier, e P, A, 
Post, G. 
Primozich, F. G. 
Pullen, E. W. 
Putnam, B. G. 


Quackenbush, R. E. 


Quinn, T. B. 
Radant, M. E. 


Raffensperger, M. J. 


Ramer, F. H., Jr 
Ramstedt, C. F. 
Rea, D. E. 

Rea, W. P., Jr. 
Redden, E. T. 
Redmond, J. G. 
Rehler, K. M. 
Rescoe, J. M. 
Rickords, T. J. 
Rieman, F. C. 
Rifkind, J. 
Robinson, J. A. 
Rogers, J. G. 
Rogers, J. M. 
Rogers, T. A. 
Romano, A. J. 
Rosenberg, R. E. 
Rosenstein, A. B. 
Ross, I. 

Ruiz, M. L. 
Russell, W. J. 
Russell, W. T. 
Saks, S. A. 
Salmi, T. W. 
Samuels, A. H. 
Sanderson, K. W. 
Sarture, C. W. 
Savant, C. J..Jr. 
Savo, T. A. 
Sawyer, H. F. 
Schalk, N. 
Schmidt, L. O. 
Schnauss, E. R. 
Schneider, R. L. 
Schneider, S. 
Schnoor, J. E. 
Schroeder, W. 
Schull, G. R. 
Schulte, R. W. 
Schultz, P. R. 
Schultz, R. T. 
Schulz, K. S. 
Scope, S. 

Scott, R. G. 
Scott, W. F. 
Scott, W. H., Jr. 
Sear, A. W. 
Sensiper, S. 
Shelley, R. G. 
Shenk, J. W. 
Shimada, G. 
Shuler, M. H. 
Shultise, Q. M. 
Shutt, S. G. 
Sibilio, R. A. 
Siegel, A Cc; 
Silva, L. M. 
Simmons, J. C. 
Simpson, B. L, 
Sink, R. L. 
Sizemore, L. E. 
Slocomb, G. M. 
Smith, B. N. 
Smith, C. C. L. 
Smith, J. D. 
Smith, L. C. 
Smith, W 
Snapp, K. 
Snyder, W. 
Sofen, I. A. 
Sohler, J. F. 
Sokol, J. L. 
Squires, W. K. 
Staudhammer, J. 
Stear, E. B. 
Steinkolk, R. B. 
Stephenson, P. L. 
Stout, T. M. 
Stow, R.A 
Sturm, T. F. 
Sundberg, R. A. 
Swanson, G. W. 
Swarthe, E. 
Takahashi, K. 


B. 
M. 
A. 


Thorenson, R. 
Tomlinson, Z. G. 
Tracey, B. P. 
Tucker, M. C, 


Van Curen, V. <2 


Vega, C.J. _ 
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Vittum, W. M. 
Vodovoz, E. 
Vulliet, P. O. 
Wachowski, H. M. 
Wakamiya, Y. 
Walker, E. S. 
Walp, R. M. 
Walsh, J. F. 
Walters, L. G. 
Wanlass, S. D. 
Watkins, E. L. 
Wedel, J. J., Jr. 
Wells, G. H 
Wennerberg, G. 
Wenters, R. L. 
Wertheim, L. A. 
Whelan, D. E. 
White, L. M. 
Whitford, R. K. 
Williams, H. M. 
Wilson, G. P. 
Wolman, L, L. 
Wong, D. S. 
Wong, E. C. 
Woo, S. H. 
Wunderlich, F. J. 
Wyatt, W. C. 
Yamada, D. A. 
Yamamoto, J. R. 
Yamamoto, T. G. 
Yoshii, H. 
Young, W. L. 
Zabusky, N. J. 
Zacharias, R. 
Ziegler, R. M. 
Zimmerman, R, L. 
Zoller, C.J. 


Phoenix 


arene Jeu 
Clark. S. 
Se ip 
Gaines, W. M. 


Ittenbach, L. i Jt. 
Levine, D. 

Mayes, TT... Laesgre 
Montgomery, E. B, 
Richman, M. A. 
Ross, J. M. 

Sabol, R. W, 
Scrafford, R. L. 
Simpson, E. J. 


Portland 


Deer, J. W. 
Doel, D 
Hinmers, K. J. 
Jenkins, R. W. 
Mendoza, D. C. 
Stone, L. N. 


Sacramento 


Akrivos, G. D. 
Cordray, R. E. 


Salt Lake City 


Hammond, S. B. 
Heiser, R. K. 
Jackson, H. L. 
Johnson, W. L. 
Linebarger, R. N. 
Lubeck, R. V. 
McCollom, Ke, Me 
Warner, B. D. 


San Diego 
Bayley, L. B., Jr. 
Biering, A. H. 
Boelens, J. 
Campbell, R. 


Fogel, L. J. 
Hallman, A. B. 
Hardy, W. G. 
Herman, J. J. 
Hodges. P. 
Holcomb, D. E., Jr. 
Kalbfell, D. C. 

La Gue, T. L. 
Loeb, M. 

Mealey, G. J. 
Michael, F. J. 
Murray, J. H. 
Norris, B. J. 
Prager, R. H. 
Scarborough, C. S. 
Schneebeck, D. A. 
Shechet, M. L. 
Stone, K. A. 
Tamura, Y. 
Waddell, B. L. 
Wade, E. 
Weisbrod, S. 


San Francisco 
Ackerlind, E. 


Allison, J. E 
Arnold, D 


Bahrs, G. S. 
Barnard, G. A., III 
Bean, D. A. 
Beatie, R. N. 
Becher, W. D. 
Bently, D. E. 
Berryhill, J. L. 
Bharucha, B. H. 
Binnall, E. P. 
Bliss) J.C; 


Boennighausen, R. A, 


Brennan, R. D. 
Bridgman, A. D., Jr. 
Brooks, H. B. 
Brooks, R. E. 
Brunetti, C. 
Bullock, J. B. 
Buntenbach, R. W. 
Bunal, T. FY 

Buss, R. R. 
Carnahan, C. W. 
Carter, J. M. 
Chesebro, E. L. 


De Tolly, GLE 8: 
Durfey, G. K. 
Ehrenburg, By Ky 
Elster, S 


Farman-Farmaian, G. 
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